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1.  SUMMARY 


1.1  Objectives 

This  two-year  project  was  directed  toward  three  major  objectives: 

•  An  assessment  of  the  event  detecticm  aitd  identification  capabilities  of  the  NORESS 
array. 

•  An  estimate  of  the  degree  to  which  these  capabilities  depend  on  the  specific  charac¬ 
teristics  of  this  site. 

•  An  evaluation  of  the  CTBT  monitoring  capabilities  using  networks  including  sta¬ 
tions  like  RSTN  and  NORESS. 


The  first  two  objectives  were  addressed  in  the  semi-annual  technical  reports  of  this  con¬ 
tract  (reports  1,  3,  and  5  in  the  following  section).  This  final  report  summarizes  these  previous 
results,  and  presents  extensive  new  work  toward  the  third  objective. 

1.2  Reports 

The  following  reports  and  papers  were  completed  under  support  by  this  contract: 

1.  Suteau-Hensoa  A.  and  T.  Bache,  Spectra  of  regional  phases  at  NORESS,  Semi- 
Anna.  Tech.  Rep.  SAIC  86II967,  Sci.  Appl.  Int  Corp.,  San  Diego,  Calif.,  1986. 

This  is  the  first  semi-annual  report  of  this  contract.  The  key  issues  addressed 
are  the  signal  characteristics  of  Pn  and  Lg  {biases  recorded  at  NORESS  (in  particu¬ 
lar,  the  repeatability  of  ^ctra  horn  blasts  in  specific  mines),  an  estimate  of  the 
average  ambient  noise  spectrum  at  NORESS,  and  an  assessment  of  event 
identification  capability  using  the  ratio  of  high  and  low  frequency  signal  energy. 

2.  Suteau-Henson,  A.,  T.  Sereno,  and  T.  Bache,  Spectral  characteristics  and  attenuation 
of  regional  phases  recorded  at  NORESS,  Proceedings  of  the  DARPA/AFGL 
Seismic  Research  Symposium,  Nantucket,  MA,  18S-190,  15-18  June  1987. 

This  short  paper  reviews  the  results  of  two  studies.  In  one  Pn  spectra  are  com¬ 
pared  for  various  events  to  evaluate  their  potential  for  event  characterization.  The 
second  smdy  presents  preliminary  results  firom  an  inversion  of  NORESS  spectra  for 
Q  and  seismic  moment. 
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3.  Sereno,  T.,  S.  Braa,  and  T.  Bache,  Regional  wave  attenuation  and  seismic  moment 
from  the  inversion  of  NORESS  spectra,  Semi-Annu.  Tech.  Rep.  SAIC  8711736,  Sci. 

Appl.  Int.  Corp.,  San  Diego,  Calif.,  AFGL-TR-87-0237,  ADA 187399,  1987. 

This  is  a  comprehensive  investigation  of  the  attenuation  of  Pn  and  Lg  spectra  j 

retarded  at  NORESS.  A  method  is  predated  for  inverting  spectra  for  both  seismic 
moment  and  apparent  attenuation.  The  method  is  applied  to  spectra  from  190 
regional  events  recorded  at  NORESS.  The  result  is  an  accurate  parameterization  of 
the  distance  and  magnitude  dependence  of  diese  spectra.  This  report  is  the  second 
semi-annual  technical  report  of  this  contract  j 

4.  Sereno,  T.,  S.  Bratt,  and  T.  Bache,  Simultaneous  inversion  of  regional  wave  spectra 
for  attenuation  and  seismic  moment  in  Scandinavia,  J.  Geophys.  Res.,  93,  2019- 
2035,  1988. 

This  paper  parallels  our  secotKl  semi-annual  technical  report  (report  3  above) 
with  the  exception  that  the  Lg  spectra  are  computed  for  a  fixed  group  velocity  win¬ 
dow  rather  than  for  a  fixed  time  length.  Our  results  for  NORESS  are  compared  to 
attenuation  estimates  used  by  Evernden  et  al.  [1986]  to  simulate  detection  capability 
in  the  Soviet  Union. 

5.  Sereno,  T.,  and  S.  Bratt,  Attenuation  atvl  detection  capability  of  regitmal  phases 
recorded  at  NORESS,  Semi-Annu.  Tech.  Rep.  SAIC  88/1598,  Sci.  Ap{4.  ItiL  Corp., 

San  Diego,  Calif.,  AFGL-TR-88-0095,  1988. 

The  main  objective  of  this  report  is  to  determine  the  relationship  between  the 
qiectral  amplitudes  used  in  our  attenuation  study  and  die  temporal  amplitudes  used 
in  signal  detection.  This  relationship  is  used  to  express  the  detection  capability  at 
NCRESS  in  terms  of  the  speural  invcrsitm  parameters  (e.g.,  Q  and  moment)  and  to 
determine  the  sensitivity  of  our  results  to  the  specific  characteristics  of  this  site.  We 
compare  our  model  to  observed  detection  statistics  and  to  the  results  of  previous 
investigations.  This  is  the  third  semi-annual  technical  report  of  this  contract. 

6.  Suteau-Henson,  A  and  T.  Bache,  Spectral  characteristics  of  regional  phases 
recorded  at  NORESS,  Bull.  Seism.  Soc.  Am..  78,  708-725,  1988. 

This  paper  parallels  the  first  semi-sumual  report  of  this  contract  (report  1  above). 

7.  Sereno,  T.,  S.  Bratt,  and  T.  Bache,  Regional  wave  detection  capability  at  NORESS, 

Proceedings  of  the  DARPA/AFGL  Seismic  Research  Symposium,  Fallbrook,  CA, 

61-66,  3-5  May  1988. 

This  short  paper  is  a  summary  of  many  of  the  results  presented  in  this  final 
report. 
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U  Suminary  of  "Estimates  of  Sdsmk  Detection  Capabiiity  in  the  Soviet  Union  Based  on 
NORESS  OhservaUons" 


13.1  Introduction 

In  this  report  we  examine  the  detection  capability  of  the  NORESS  array  and  its  implica¬ 
tions  for  estimating  the  detection  threshold  of  a  network  of  seismic  stations  in  and  around  the 
Soviet  Union.  NORESS  was  designed  as  a  prototype  array  for  regiorud  monitoring,  and  it  is 
located  within  regional  distance  to  portions  of  the  western  USSR.  Therefore,  the  NORESS 
detection  capability  provides  a  reasonable  basis  for  normalizing  simulations  of  the  detection 
capability  of  hypothetical  networks  in  the  vicinity  of  the  Soviet  Union.  However,  NORESS 
may  not  represent  the  actual  attenuation  and  noise  for  the  stations  in  our  hypothetical  network, 
so  we  also  determine  the  sensitivity  of  the  results  to  changes  in  the  frequency-dependent  signal 
and  noise  characteristics.  The  primary  issues  ittvolved  in  seismic  monitoring  ate  detection, 
location,  and  event  identihcation,  and  these  are  discussed  below. 

•  Detection  Capability.  The  primary  focus  of  this  report  is  on  the  detection  capability  of 
hypothetical  rretworks  in  the  Soviet  Union. 

•  Location.  Signal  detection  is  only  meaningful  if  there  is  adequate  information  to  locate 
the  source.  Therefore,  our  simulations  of  detection  capability  require  detection  of  at  least 
three  phases  (arrival  times  and  azimuths)  by  the  hypothetical  network  (constraining  lati¬ 
tude,  longitude,  and  depth).  Bratt  et  al.  (1987]  show  that  for  a  rretwork  irKluding  20 
NORESS'type  arrays  internal  to  the  Soviet  Union  and  13  external,  the  3-phase  location 
uncertainty  for  events  near  the  detection  threshold  is  about  10-15  km.  This  uncertainty 
can  be  reduced  with  more  accurate  travel-time  tables  and/or  master-event  location  tech¬ 
niques,  as  discussed  in  that  report.  Also.  Bratt  et  al.  [1987]  analyzed  the  sensitivity  of 
the  location  uncertainty  to  network  configuration.  We  present  no  further  work  on  these 
issues  here. 

•  Event  Identification.  Numerous  techniques  have  been  proposed  to  discriminate  between 
explosions  and  earthquakes  at  regional  distances  [for  review,  see  Pomeroy  et  al ,  1982], 
It  a(^ars  that  regiorud  event  identification  requires  the  ability  to  distinguish  relatively 
subtle  source  effects  within  a  cloud  of  complex  propagation  effects.  For  this  reason,  it  is 
important  to  accurately  calibrate  propagation  effects  for  specific  regions  of  interest. 

In  this  report  we  concentrate  on  detection  capability  because  it  is  the  most  fundamental 
aspect  of  mcmitoiing  capability.  Lowering  the  detection  threshold  gives  roughly  proportional 
improvements  to  the  location  and  identification  capabilities. 

13.2  Approach 

Our  approach  to  estimating  the  detection  capability  of  a  hypothetical  network  in  the 
Soviet  Union  involves  the  following  steps; 
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1.  We  drennine  an  accurate  parameterization  of  regional  Pn  and  Lg  spectra  recorded  at 
N^-RESS  in  terms  of  evettt  magnitude  and  apparent  attenuation.  This  is  done  by  invert¬ 
ing  spectra  rimn  186  regional  events  for  seismic  moment  and  apparent  attenuation.  With 
this  parameterization,  we  can  predict  the  spectra  for  a  chosen  source  and  range  and  be 
conlidem  that  our  results  are  conect  vrithin  the  well-defined  uncertainty  of  our  parame¬ 
ters. 

2.  We  determine  the  relationship  between  the  spectral  amplitudes  aixl  the  temporal  ampli¬ 
tudes  used  in  sigiud  detection.  This  allows  us  to  use  our  spectral  parameterization  to  esti¬ 
mate  (temporal)  detection  ctqiability. 

3.  We  estimate  NORESS  detection  cryjability  with  the  spectral  parameterizatioa  Since  this 
detection  capability  is  a  function  of  the  attenuation  observed  at  NORESS  and  other  terms 
specific  to  the  NORESS  array  configuration  artd  processing,  we  can  estimate  the  extent  to 
which  our  results  depend  on  characteristics  specific  to  NORESS.  Our  NORESS  detection 
capability  model  is  validated  with  observed  detection  statistics. 

4.  We  simulate  the  detection  capability  of  a  hypothetical  network  of  NORESS-capability 
arrays  in  the  Soviet  Union.  Of  course,  this  assumes  that  NORESS  capability  represents 
the  signal  and  noise  environment  in  die  Soviet  Union.  Therefore,  we  estimate  the  sensi¬ 
tivity  of  these  simulations  to  changes  in  the  assumptions  about  the  attenuation  and  noise 
in  the  Soviet  Uruon. 

13.3  Results 

We  have  determined  a  simple  and  accurate  parameterization  of  the  magnitude  and  dis¬ 
tance  dependence  of  regional  Pn  and  Lg  spectra  recorded  at  NORESS.  This  parameterization 
is  based  on  the  inversion  of  spectra  from  186  events  with  magnitudes  between  1.1  and  4.8  and 
epicentral  distances  between  2(X)  and  1400  km.  Ingure  1.1  displays  the  predicted  Pn  and  Lg 
spectra  at  500,  8(X),  and  1000  km  for  a  magnitude  3.0  event,  based  on  this  parameterization. 
The  bottom  curve  is  the  average  ambient  noise  spectrum  at  NORESS  [Suteeut-Henson  and 
Bache,  1988].  The  predicted  Pn  spectrum  is  nearly  parallel  to  the  noise  curve  at  a  range  of 
about  500  km,  a  result  that  is  consistent  with  a  NORSAR  study  that  examined  data  from  the 
high-frequetKy  element  of  the  NORESS  array  [Ringdal  et  al.,  1986],  At  longer  ranges,  the 
predicted  frequetKy  of  the  maximum  signal-to-noise  ratio  is  considerably  less  than  15  Hz. 
Evernden  et  al.  [1986]  propose  that  detection  capability  would  be  enhanced  at  high  frequency 
(>  20  Hz).  Our  results  for  Scandinavia  do  not  support  this  contention  except  at  ranges  less 
than  about  400-500  km. 

Figure  1.1  is  an  accurate  represeruation  of  Lg  spectra,  but  cannot  be  used  to  estimate 
signal-to-noise  (the  ambient  noise  curve  is  plotted  with  the  Lg  spectra  in  Figure  1.1b).  The 
"noise”  prior  to  the  Eg  arrival  is  actually  the  ambient  noise  plus  P  and  Sn  coda,  and  is  a  com¬ 
plicated  function  of  range  and  source.  However,  we  can  note  that  the  Lg  spearum  reaches  the 
ambient  noise  level  at  a  much  lower  frequency  than  the  Pn  spectrum,  while  at  long  periods  the 
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Lg  ampiitiirif;  exceeds  that  of  Pn  by  as  much  as  a  factor  of  10.  This  is  consistent  with  the 
observation  that  regional  seismograms  are  characterized  by  Lg  being  the  largest-amplitude 
phase  and  Pn  having  a  higher  dominant  frequency. 

To  pstimatp.  detection  capability  from  these  spectral  results,  it  is  necessary  to  determine 
the  relationship  between  temporal  and  spectral  amplitudes.  Signal  detection  at  NORESS  is 
achieved  when  the  short  term  average  (STA)  divided  by  the  long  term  average  (LTA)  exceeds 
a  predetermined  threshold.  These  amplitudes  are  measured  on  filtered  beams.  Therefore,  it  is 
necessary  to  determine  the  relationship  between  single-channel  spectral  amplitudes  and  these 
time-domain  amplitudes  on  filtered  beams.  For  Pn  we  find  that  the  temporal  STA/LTA  can  be 
expressed  as  a  product  of  the  single  channel  spectral  SNR  (for  which  we  have  parameterization 
in  terms  of  Q  and  magnitude)  and  terms  specific  to  the  array  configuration  and  beamforming. 
This  factorization  allows  us  to  determine  the  extent  to  which  our  NORESS  Pn  detection  capa¬ 
bility  estimates  depend  upon  the  specific  characteristics  of  this  site. 

The  Lg  spectrum  depends  on  wirulow  length  and  the  prc-Lg  noise  includes  the  coda  of 
previously  arriving  phases.  Therefore,  a  similar  relationship  can  irot  be  found  for  Lg.  We  sim¬ 
ply  note  that  the  Lg  spectral  amplitudes  computed  for  5-s  windows  starting  just  prior  to  the 
detection  time  ate  proportional  to  the  STA  on  the  beam  with  the  maximum  SNR.  However, 
Lg  Q  estimates  based  on  fixed  5-s  windows  are  trot  available  for  other  regions,  so  our  results 
for  Lg  are  not  easily  extrapolated.  The  prc-Lg  LTA  is  more  complex.  We  parameterized  it  as 
a  function  of  frequency,  distance,  beam  type  (coherent  or  incoherent),  and  magnitude,  and  the 
latter  is  most  important  We  fitrd  that  the  log  LTA  iircteases  proportional  to  0.8  Since 
both  Lg  STA  and  LTA  increase  with  fire  detectability  of  Lg  is  much  less  dependent  on 
magnitude  than  the  detectability  of  Pn. 

We  estimate  NORESS  detection  capability  using  our  spectral  inversion  results  together 
with  an  empirical  parameterization  of  the  temporal  SNR.  The  latter  is  called  the  "detectabil¬ 
ity,"  and  it  gives  the  probability  of  detection  of  a  given  magnitude  event  as  a  function  of  epi- 
central  distance.  The  detectability  is  computed  directly  from  detector  parameters  (e.g.,  the 
STA,  LTA,  and  pre-set  detector  thresholds).  For  extrapolating  these  results  to  other  areas 
(with  estimates  for  the  attenuation)  we  use  our  spectral  inversion  results  and  the  relationship 
between  temporal  and  spectral  amplitudes.  To  validate  this  procedure,  we  compare  "predicted” 
(based  on  our  inversion  results)  and  "observed”  Pn  detectability  at  NORESS.  The  dashed 
curves  in  Figure  1.2  plot  our  estimates  of  the  90%  NORESS  thresholds  for  detecting  Pn 
and  for  detecting  Pn  or  Lg,  based  on  the  "detectability."  The  solid  curve  is  the  90%  thres¬ 
hold  for  detecting  Pn  based  on  the  inversion  results.  Also  indicated  in  Figure  1.2  are  the 
results  of  Ringdal  [1986]  who  estimated  the  90%  threshold  for  detecting  P  phases  between 
700  and  14(X)  km  to  be  2.7  and  for  detecting  P  or  secondary  phases  in  the  same  distance  range 
to  be  2.5.  The  average  distance  of  the  events  used  in  that  study  is  about  960  km.  The  three 
independent  studies  produce  consistent  estimates  of  the  regional  wave  detection  capalxlity  of 
the  NORESS  array,  which  lends  considerable  support  to  our  spectral  parameterizatioa 

Our  confidetKe  in  predictions  of  detection  capability  for  other  regions  depends  on  how 
well  we  know  the  regional  attenuation  and  rwise.  Simulating  the  detection  c^^>ability  of 
hypothetical  networks  within  the  Soviet  Union  is  difficult  because  accurate  attenuation  esti¬ 
mates  and  noise  spectra  are  not  available  for  that  region.  Therefore,  we  use  a  range  of  models 
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Figure  1.2.  Estimates  of  90%  detection  thresholds  at  NORESS  as  a  function  of  epicentral 
distance.  The  dashed  curves  are  based  on  a  parameterization  of  observed  temporal  amplitudes. 
Curves  for  detecting  Pn  and  for  detecting  either  Pn  or  Lg  are  plotted.  The  solid  curve  is  based 
on  the  results  of  our  inversion  of  Pn  spectra.  The  horizontal  bars  are  detection  thresholds  at 
NORESS  as  determined  by  comparing  NORESS  detections  to  bulletins  {Hoduced  by  local 
seismic  networks  [Ringdal,  1986].  Events  for  this  study  were  at  ranges  between  700  and  1400 
and  the  average  distance  was  about  960  km  (X). 
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in  our  netwoik  capaMity  simulations  to  determine  the  sensitivity  of  our  results  to  attenuation 
and  noise  characteristics.  As  an  example.  Figure  1.3  shows  a  SNAP/DX  simulation  [Bratt  et 
al.,  1987]  of  the  90%  threshold  for  detecting  three  Pn  phases  from  a  netwoik  composed  of 
NORESS-quality  stations;  13  outside  and  20  inside  the  USSR.  In  this  case  the  NORESS 
detectability,  including  both  attenuatkm  and  noise,  was  used  to  represent  Pn  propagation  and 
detection  at  all  stations.  That  is.  Figure  1.3  represents  the  detection  capability  if  signal  and 
noise  are  identical  to  that  observed  at  NORESS  throughout  the  Soviet  Unioa  In  teims  of 
NORESS  Ml,  the  threshold  is  between  2.4  and  2.7  in  die  Soviet  Union.  A  S0%  increase  in  Pn 
Q  reduces  the  Mi  threshold  to  between  2.3  and  2.6,  and  a  50%  decrease  in  Q  increases  the  Mi 
threshold  to  between  2.7  and  2.9.  Including  Lg  reduces  the  Mi  threshold  by  about  0.2  to  0.3. 
To  improve  confidence  in  these  simulations  of  network  capability,  it  will  be  necessary  to  obtain 
accurate  estimates  of  the  frequency-dependence  of  atteiuation  and  noise  in  the  Soviet  Union. 

13.4  Outline  of  report 

This  report  includes  five  technical  sections.  Section  2  reviews  die  method  and  results  of 
our  generalized  inversion  of  regional  wave  spectra.  This  is  a  summary  of  material  in  previous 
reports  [Sereno  et  al.,  1987;  Sereno  and  Bratt,  1988].  Section  3  describes  the  relationship 
between  temporal  and  qiectral  amplitudes  and  presents  our  estimates  of  the  detection  capability 
at  NORESS.  This  is  a  summary  of  the  material  presented  in  our  last  semi-annual  report 
[Sereno  and  Bratt,  1988]. 

Section  4  is  new  work  presented  for  die  first  time  on  detection  capability  in  the  Soviet 
Union.  Section  4.1  discusses  die  procedures  and  uncertainties  involved  in  extrapolating  the 
NORESS  results  to  simulate  the  detection  capability  of  a  network  of  seismic  stations  in  the 
Soviet  Union.  Section  4.2  presents  detection  capability  simulations  for  a  variety  of  assump¬ 
tions  about  the  attenuation  and  nrfise  in  the  Soviet  Union.  Section  5  summarizes  our  major 
conclusions  developed  during  the  course  of  this  contract,  widi  emidiasis  on  the  new  results  in 
Section  4. 
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Figure  1.3.  COTlours  of  the  90%  threshold  for  detecting  3  Pn  phases  for  a  network  of 
NORESS-type  arrays;  13  external  to  the  Soviet  Union  and  20  internal.  It  is  assumed  that 
attenuation  and  noise  are  identical  to  that  observed  at  NORESS. 


2.  SPECTRAL  INVERSION 


An  important  oonsideraiion  in  designing  a  seismic  netwoilc  intended  to  monitor  com[^- 
ance  with  a  miclear  test  ban  treaty  is  the  spectral  content  of  regional  wave  signals.  For  exam¬ 
ple,  Evernden  et  al.  [1986]  aigue  that  efficient  Ingh-fiequency  (>  20  Hz)  Pn  propagation  can 
be  exploited  to  reduce  detection  thresholds  in  the  Soviet  Union  and  to  aid  in  event 
identification.  Of  course,  this  im(Aies  that  high-frequency  sensors  should  be  an  integral  part  of 
the  seismic  network.  However,  fiieir  argument  was  based  primarily  on  1-3  Hz  amplitude  data 
from  paths  in  the  eastern  United  States,  and  the  only  high-frequency  data  they  use  to  corro¬ 
borate  their  result  is  a  single  Pn  qiectnim  recorded  at  a  range  of  only  190  km.  In  this  section 
we  address  diis  issue  with  a  large  data  set  This  is  a  summary  of  the  method  and  results  of 
our  previous  study  of  frequency-dependent  regional  wave  propagation  for  paths  to  NORESS 
[Sereno  et  al.,  1987,  Sereno  et  al..  1988].  In  these  studies,  we  inverted  spectra  from  186 
regional  events  covering  distances  from  200  to  1400  km  aitd  magnitudes  between  1.1  and  4.8 
for  both  seismic  moment  and  {4)parent  attenuation.  Pn  spectra  were  inverted  between  1-lS  Hz 
and  Lg  ^ctra  betweoi  1-7  Hz.  The  final  result  is  an  accurate  parameterization  of  the  range 
and  magnitude  dependence  of  these  spectra.  We  note  that  our  results  for  NORESS  are  incon¬ 
sistent  with  the  attenuatkm  model  used  by  Evernden  et  al.  [1986],  and  they  do  not  support 
their  contention  of  enhanced  detection  capability  at  high  frequency,  except  at  rsmges  less  than 
400-500  km. 

2.1  Method 

This  sectiem  is  a  brief  summary  of  the  method  we  used  to  estimate  attenuation  and  source 
parameters.  A  detailed  description  is  given  by  Sereno  et  al.  [1988].  The  analysis  assumes  a 
simfde  source  scaling  model  and  that  all  observed  spectra  can  be  fit  by  a  single  hequency- 
dependent  Q  model  We  parameterize  the  instrument-correaed  amplitude  spectrum  of  a 
seismic  signal  as 


A(f,r)  -  S(f)  G(r,ro)  exp 


-1^ 

m 


(2.1) 


where  A(f/)  is  the  observed  disi^acement  spectrum  at  range  r  and  frequetKy  /,  Slf)  is  the 
source  spectrum,  G{r,r^  is  geometric  spreading,  and  the  last  term  is  the  effective  attenuation 
for  travel  time  t.  The  effective  attenuation  includes  contributions  from  both  anelasticity  and 
scattering. 

Source  spectra.  We  adopt  a  simi^ified  Mueller  and  Murphy  [1971]  explosion  source 
function  characterized  by  f~^  decay  beyond  a  conrer  frequency  that  scales  inversely  with  the 
cube  root  of  the  long-period  level.  While  source  assumptions  have  had  a  controlling  influence 
in  some  studies  of  attmuation,  they  are  not  crucial  for  our  analysis  since  most  of  the  events 
had  magnitudes  less  than  3.0  and  apparent  corner  ffequencies  beyond  the  frequency  band 
inverted.  The  source  parameters  estimated  by  the  inversion  are  the  long-period  level,  Sq,  for 
each  evoit  and  the  corner  frequency  scaling  parameter.  The  explosion  moments  are  estimated 
ftom  the  long-period  levels  derived  from  Pn  [Stevens  and  Day,  1985].  Since  our  observations 
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are  fiom  a  single  station  and  the  focal  medianisms  for  the  small  earthquakes  in  our  data  set  are 
unknown,  we  do  not  know  the  radiatitm  pattern  and  are  therefore  unable  to  estimate  seismic 
moment  fiom  the  earthquake  source  level  derived  from  Pn.  However,  since  Lg  samples  a  large 
fraction  of  the  focal  qrhere,  the  earthquake  moments  can  be  estimated  from  the  long-period 
source  levels  derived  fiom  Lg  [Street  et  al..  1975].  The  relationships  we  use  to  estimate 
seismic  moment  are 

(2.2) 

=4jtp^^(Lg)  (2.3) 

where  p,  and  a,  are  near-surface  density  and  compressional  velocity,  and  ^  average 
crustal  density  and  shear  wave  velocity,  and  exp  and  eq  denote  explosions  and  earthquakes, 
respectively. 

The  amount  of  Lg  energy  excited  by  an  explosion  is  depth-dependent  and  complicated  by 
near-source  wave  conversions  [Bennett  et  al.,  1987].  Explosions  generate  Lg  energy  primarily 
through  P-SV  mode  conversions  and  scattering,  while  earthquakes  produce  much  more  direct 
shear  wave  energy.  Therefore,  for  a  given  source  moment  and  focal  depth,  eardtquakes  are 
expected  to  generate  larger  Lg  amplitudes  than  explosions.  We  will  express  the  long-period  Lg 
source  level  for  an  explosion  as  an  unknown  fraction  of  the  long-period  level  for  an  earthquake 
of  equal  moment  and  depth.  That  is, 

S§^(Lg)  =  KSSf«(Lg)  (2.4) 

where  K  is  an  unknown  constant,  presumably  less  than  one.  We  estimate  this  constant  using 
(2.2),  (2.3),  and  (2.4),  together  with  the  observed  ratio  of  S^(Lg)  to  S^(Pn). 

Geometric  spreading.  The  inversion  also  requires  that  we  assume  the  geometric  spread¬ 
ing  fimction.  FoUowing  Herrmann  and  Kijko  [1983],  we  express  the  spreading  function  in  the 
frequency  domain  as 

G{r,r^  -  (1/^)  for  r  ^  r^ 

(ro/r)*"  for  r  ^  Tq 

where  ro  is  a  transition  distance  from  spherical  spreading  to  spreading  rate  m.  The  Lg  phase 
consists  of  higher-mode  surface  waves  which  are  accurately  described  at  long  ranges  by 
cylindrical  spreading  (m  s  1/2),  provided  the  window  length  is  sufficient  to  encompass  the 
entire  dispersed  wave  train.  By  comparing  the  long-period  amplitude  spectrum  of  Lg  to 
moments  calculated  from  long-period  surface  waves.  Street  et  al.  [1975]  empirically  deter¬ 
mined  cq  =  100  km,  or  roughly  twice  the  crustal  thickness.  Measuring  the  decay  rate  of  syn¬ 
thetic  Lg  phases  computed  for  an  elastic  medium,  Herrmann  and  Kijko  [1983]  verified  that  Lg 
frequency  domain  spreading  was  accurately  described  as  cylindrical  and  substantiated  the 
empirical  result  of  Street  et  al.  [1975]  for  tq.  Therefore,  we  adopt  (2.5)  with  tq  =  100  km  and 
m  =  1/2  to  approximate  Lg  geometric  spreadirig. 
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Less  woik  has  been  done  on  die  spfeading  rate  of  Pn.  Because  its  energy  density  is 
more  localized  about  a  single  ray  path,  Pn  geometric  qifeading  is  more  sensitive  to  velocity 
gradients  in  the  upper  raarttle.  Numerical  studies  of  Pn  indicate  that  for  typical  upper  mantle 
structures  the  spreading  rate  lies  between  r~'  and  r~^  [Langston,  1982;  T.  C.  Wallace,  persrmal 
cmiimunication,  1987].  We  find  that  an  importam  constraint  on  Pn  spreading  is  the  con¬ 
sistency  of  the  derived  source  parameters  from  the  separate  Lg  and  Pn  inversions.  This  cri- 
terum  supports  a  choice  of  for  the  Pn  qireading  rate.  Therefore  we  use  (2.5)  with  =  1 
km  and  m  =  1.3  to  describe  Pn  spreading.  We  experimented  with  odier  spreading  rates  and 
found  only  minor  differences  inQ.  A  chan^  in  the  assumed  sfneading  rate  simply  trades-off 
with  derived  moment  since  the  observations  are  fimn  a  arigle  station. 

Effective  attenuation.  We  characterize  the  range-dependeitt  decay  of  the  seismic  spectrum 
in  terms  of  a  power  law  frequency  dependence  of  Q.  That  is. 

Qif)  *  Qo/”  (2.6) 

where  (2o  H  are  parameters  of  die  inversitm.  No  attempt  is  made  to  distinguish  intrinsic 
absorption  from  scattering.  In  this  form,  our  results  are  easUy  compared  to  those  for  other 
geogiaidiic  regions.  We  have  not  accounted  for  any  azimuthal  variations  in  Q,  but  have  com- 
Inned  dtua  from  all  azimuths  into  an  inversion  for  a  sirigle.  average  Q  model.  Examination  of 
the  fit  of  the  model  to  the  data  demonstrates  the  validity  of  this  procedure. 

Inversion.  The  input  data  for  die  inversim  are  the  logarithms  of  the  observed  displace¬ 
ment  spectra  corrected  for  die  assumed  geometric  spreading.  These  data  are  inverted  for 
^iparent  attenuatkm.  source  moment,  and  the  constant  relating  comer  frequency  and  moment 
Adrqptmg  standard  methods  for  solving  non-linear  inverse  pcoUems,  we  linearize  the  system  of 
equations  governing  the  relaricmship  between  the  data  and  model  parameters.  We  assume  a 
starting  modd,  compute  theoretied  data,  subtract  it  fiom  the  observed  data,  and  solve  itera¬ 
tively  fm*  the  model  peiturtMtions  that  minimize  the  data  residual  in  die  least  squares  sense.  In 
practice,  we  have  found  it  necessary  to  include  damping  to  stablize  the  solution.  That  is,  we 
minimize  a  weighted  sum  of  the  data  residuals  and  the  model  perturbation  norm.  The  explicit 
problem  formulation,  matrices,  and  partial  derivatives  are  given  by  Sereno  et  al.  [1987]. 

2.2  Data 

The  data  used  in  this  study  consist  of  staUe,  array-averaged  spectra  for  186  regional 
events  recorded  by  the  small  aperture  NC^IESS  seismic  array  in  Norway.  Event  magnitudes 
range  between  1.1  and  4.8  and  epicentral  distances  are  between  2(X)  and  14(X)  km.  The  array 
includes  25  short-period  instruments  in  concentric  rings  witti  a  maximum  diameter  of  3  km. 
The  data  are  digitally  recorded  at  40  samples  per  second.  The  NORESS  short-period  instru¬ 
ment  response  is  qipioximately  flat  to  velocity  between  1  and  10  Hz. 

Signal  processing.  The  Pn  spectra  were  calculated  as  part  of  an  automated  seismic  array 
processing  program  (SAIAP)  developed  as  an  extension  of  the  RONAPP  program  used  at  NOR- 
SAR  [Mykkeltveit  and  Bungum,  1984].  The  program  computes  spectra  for  each  automatically 
detected  signal.  The  spectral  estimation  technique  is  that  proposed  by  Bache  et  al.  [1985].  A 
10%  cosine-squared  taper  is  tqiplied  to  a  5-s  window  starling  0.3  s  before  the  onset  time  of  the 


arrival  cm  the  veitical  OMnponem.  The  time  series  is  padded  with  zeros  to  512  samples  and 
fast  Fourier  transformed.  The  same  pfoceduie  is  aj^lied  to  a  noise  sample  taken  prior  to  the 
first  P  detection.  The  squared  noise  amplitude  spectrum  (power)  is  subtracted  from  the 
squared  signal  spectrum  (eneigy  density).  The  resulting  noise-corrected  signal  spectra  are  aver¬ 
aged  across  the  array  and  corrected  for  the  instniment  response.  Bache  et  al.  [1985]  show  that 
if  the  noise  is  random,  stationary,  and  uncorrelated  with  the  signal,  the  signal  spectrum  esti¬ 
mate  obtained  with  this  method  converges  to  the  tnte  signal  spectmm  as  the  number  of  ele- 
meats  increases.  Array  averaging  also  has  the  desirable  effect  of  suppressing  uncorrelated 
local  site  effects.  We  experimented  with  Pn  window  lengths  of  5,  10,  and  15  s  and  found  that 
the  spectra  are  insensitive  to  that  parameter. 

The  inversion  results  presented  here  were  obtained  using  Lg  spectra  computed  for  a  fixed 
group  velocity  window  of  3.6-3.0  km/s.  The  spectra  were  computed  for  each  array  element, 
corrected  for  the  ambient  (pre-/*n)  noise,  and  array-averaged  using  the  method  of  Bache  et  al. 
[1985].  Our  low  group  velocity  cutoff  was  chosen  as  3.0  km/s  because  our  events  typically 
produced  Lg  phases  with  low  signal-to-noise  for  group  velocities  less  than  this.  Fixed  group 
velocity  windows  ensure  that  the  same  modes  contribute  ro  the  Lg  spectnim  at  aU  epicentral 
distances.  The  Q  estimated  from  these  spectra  is  a  measure  of  the  average  absorption  of  shear 
waves  in  the  crust  [e.g.,  Campillo  et  al.,  1985].  However,  the  energy  distributed  over  such  a 
long  time  window  (56  s  at  KXX)  km)  is  not  an  appropriate  measure  for  detection  capability 
assessment  We  find  the  ^)ectra  computed  for  short  fixed  time  lengths  are  more  appropriate 
for  that  application.  Sereno  and  Bratt  [1988]  compare  the  inversion  results  using  Lg  spectra 
computed  with  fixed  group  velocity  and  fixed  time  length  windows. 

Events.  Figure  2. 1  {dots  the  epiceriters  of  all  of  the  events  used  in  the  inversion.  A  com¬ 
plete  list  of  locations,  origin  times,  and  magnitudes  are  given  by  Sereno  and  Bratt  [1988].  The 
location  and  origin  times  are  from  a  local  bulletin  published  by  the  University  of  Bergen  or  the 
University  of  Helrinki,  when  available,  or  from  the  Preliminary  Determination  of  Epicenters 
(PDE)  bulletin.  Events  for  which  an  indepeiKlent  network  solution  is  not  available  are 
assigned  SAIAP  locations.  The  are  based  on  the  Lg  amplitude  computed  by  RONAPP. 
They  differ  slightly  from  the  RONAPP  NORESS  bulletin  by  being  distance-corrected 

to  the  event  location  computed  by  one  of  the  independent  networks. 

Mining  explosions  are  rqxrrted  in  the  Helsinki  bulletin  and  by  the  University  of  Bergen. 
Events  that  are  not  reported  explosions,  but  have  locations  within  50  km  of  a  Iroown  mine  are 
considered  of  unknown  source  type.  Other  events  are  presumed  to  be  earthquakes.  The  186 
events  include  107  explosions,  63  presumed  earthquakes,  and  16  events  of  unlroown  source 
type.  In  some  cases  only  one  phase  was  included  because  the  other  was  not  detected  or 
because  it  had  a  low  signal-to-noise  ratio  over  the  frequency  band  used  in  the  inversion.  The 
Pn  inversion  included  152  events,  with  83  explorions,  56  presumed  earthquakes,  and  13  events 
of  unktrown  soutx%  type.  The  Lg  inversion  used  160  events,  including  92  explosions,  53 
presumed  earthquakes,  and  15  unknown.  Of  the  186  events,  126  were  used  for  both  the  Pn 
and  the  Lg  analyses. 
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23.1  Seismic  moment  and  Q 

For  a  fixed  set  of  source  and  spreading  assumptions,  the  separate  Pn  and  Lg  inversions 
each  define  a  range  of  models  consistent  with  the  observations.  For  example,  the  Pn  inversion 
produced  models  with  Q  between  200  and  SOO that  differed  in  data  variance  by  less 
than  5%.  Similarly,  Lg  Q  models  between  420 and  570  produce  data  variances  that 
differ  by  less  than  2%.  To  resolve  these  parameter  trade-offs,  we  impose  the  constraint  that 
the  separate  Pn  and  Lg  inversions  give  an  internally  consistent  set  of  source  parameters.  In 
particular,  we  require  that  the  ratio  of  the  Lg  and  Pn  source  levels  be  independent  of  epicentral 
distance.  This  reduces  the  range  of  acceptable  models  to  a  set  of  model  pairs  consistent  with 
both  sets  of  observations.  Table  2.1  lists  examples  of  these  model  pairs.  Note  that  none  of  Lg 
models  are  consistent  with  Pn  Qq  greater  than  350. 


Table  2.1.  Q  Models  With  Internally  Consistent  Source  Parameters 


Qi^ 

Qpit 

420 /®^ 

225 

500/®^ 

280/®^^ 

560/®^ 

325  /®^ 

There  is  no  obvious  basis  for  choosing  among  these  model  pairs.  However,  Sereno  et  al. 
[1988]  note  that  a  trade-off  exists  betweoi  data  variance  pmJ  variance  in  the  moment- 
magiutude  relationship.  Based  on  this  trade-off.  they  selected  Q^g(/)  =  560 as  their  "pre- 
ferred"  model.  Applying  the  source  consisteix;y  constraint,  tiiis  supports  a  Pn  Q  model  with 
Qo  -  325  and  q  =  0.48.  Although  this  trade-off  analysis  is  subjective,  it  is  only  used  to 
choose  among  a  relatively  small  range  of  Ql,  models. 

The  earthquake  moments  are  estimated  from  the  source  levels  derived  fn>m  Lg  using 
(2.3)  with  Pc  s  2.7  gm/cm^  and  =  3J  km/s.  The  explosion  rnmnents  are  derived  from 
S^fPn)  using  (2.2)  with  p,  »  2.5  gm/cm^  and  =  5.0  km/s.  Figure  22  displays  the  derived 
source  moments  as  a  function  of  NORESS  The  solid  lines  indicate  the  least  squares  linear 
fits  given  by 


log  =  1.08  Ml  17.6  (2.7) 

log  M^  =  1.03  Ml  +  17.1  (2.8) 

However,  because  of  an  uneven  distribution  of  magnitudes  in  our  explosion  data  set,  the  least 
squares  fit  does  not  adequately  represent  the  observed  Mi  dependence  of  log  Mq.  The  fit 
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Pn  Inverted  Moment  vs  M|^  -  Explosions 


Lg  Inverted  Moment  vs  M|^  *  Earthquakes 


FiguR  2.2.  Invened  oNiBem  veisus  NORESS  local  mapiitude  for  (a)  explosiaas  and  (b) 
eaithqaakes.  The  ejqdocioa  momentt  were  derived  ftoo  Pn  and  the  eafthquake  momenii  were 
derived  fh»  Lg.  The  leaat-sqnares  linear  fit  of  log  A#o  to  Af^  is  {dotted  as  the  solid  line  in 
each  Ague.  The  dashed  line  in  Figure  2.2a  |dou  the  linear  trend  for  eiqdosions  fitted  by  eye. 
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underestimates  the  moments  of  large  events  and  overestimates  the  moments  of  small  events. 
Therefore,  we  simply  fit  the  linear  trend  by  eye.  The  result  is  the  dashed  line  in  Figure  2.2 
which  is  expressed  as 


log  A/?**  =  1.35  Ml,  +  16.90  (2.9) 

The  coefficient  of  Mi  is  higher  than  that  expected  if  Mi  is  approximately  equal  to  mi,.  The 
explosion  moments  derived  from  Lg  spectra  display  a  similar  rate  of  increase  with  Mi,  while 
the  earthquake  log  moments  display  much  less  scatter  and  increase  at  a  rate  close  to  1.0  Af^ 
(Figure  2.2b).  Therefore,  it  is  not  likely  that  inaccurate  attenuation  parameters  are  responsible 
for  the  high  Mi  coefficient  observed  for  explosions.  Since  Mi  is  based  on  Lg  amplitude,  it  is 
possible  that  the  high  Mi  coefficient  is  due  to  complexity  in  Lg  excitation  by  explosions. 

Sereno  et  al.  [1988]  determine  the  sensitivity  of  the  inversion  results  to  assumptions 
about  the  Pn  geometric  spreading  rate.  Spreading  rates  of  r~^,  and  were  explicitly 
investigated,  and  for  each  there  is  a  Q  model  that  can  adequately  reproduce  the  Pn  spectra. 
However,  since  our  parameterization  of  geometric  qireading  is  fiequency-independent  and  our 
dau  is  from  a  single  station,  changes  in  the  assumed  spreading  rate  trade-off  directly  with  fire 
source  moment,  and  have  very  little  affect  on  the  range  of  Q  models  produced  by  the  Pn  inver¬ 
sion.  However,  an  important  constraint  on  the  Pn  spreading  rate  is  the  relative  excitation  of  Lg 
for  explosions  and  earthquakes  of  equal  mometu.  That  is.  since  ic  in  (2.4)  is  estimated  using 
the  explosion  moments  derived  fttm  Pn.  the  most  dramatic  effect  of  changing  the  Pn  spreading 
rate  is  the  implied  change  in  the  relative  Lg  exciution  of  explosions  and  earthquakes.  We  find 
that  for  earthquakes  and  exfdosions  of  equal  moment,  ic  »  0.27  for  our  "preferred"  attenuation 
model  This  means  that  the  average  Lg  earthquake  excitation  is  appmximately  four  times  the 
average  Lg  explosion  excitation  for  sources  of  equal  moment  This  is  consistent  with  results  of 
previous  studies  that  compared  Lg  amfditudes  from  approximately  colocated  earthquakes  and 
explosions  [e.g.,  WtUis,  1963;  Pomeroy,  1977;  NuttU,  1981].  Assumed  Pn  spreading  rates 
much  differem  from  give  earthquake  to  explosion  Lg  excitation  ratios  which  are  incon- 
sistem  with  these  empirical  observations  [see  Sereno  et  al.,  1988]. 

Table  2.2  surrunarizes  the  inversion  results  for  our  "preferred"  model  The  only  parame¬ 
ter  not  listed  is  comer  frequency.  Comer  finequency  was  assumed  to  scale  inversely  with  the 
cube-ro(X  of  the  moment  Our  results  indicate  that  an  Mi  =  3.0  event  has  an  approximate 
comer  fiequency  of  10  Hz,  although  our  data  do  not  accurately  resolve  that  parameter.  How¬ 
ever,  we  note  that  our  results  are  within  a  few  hertz  of  those  from  near-field  estimates  for 
events  with  the  same  local  magnitude  [Bungum  et  al.,  1982;  Mueller  and  Cranswick,  1985]. 
IHgure  2.3  shows  some  examples  of  the  fit  of  our  model  to  the  observed  spectra.  The  spectra 
are  plotted  over  the  bandwidth  inverted;  1-15  Hz  for  Pn  and  1-7  Hz  for  Lg.  These  examples 
are  typical  of  the  quality  of  the  fits  that  were  obtained  for  most  of  the  events.  A  complete 
catalog  comparing  observed  and  predicted  ^ctra  is  given  by  Sereno  and  Bran  [1988]. 
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Figure  2.3.  Selected  examplee  of  the  fit  of  Pn  and  Lg  apectia  bwed  on  the  invenion  resulti  to 
obaenred  spectre.  The  Pn  spectre  are  plooed  on  the  left  and  the  Lg  spectre  on  the  right  The 
spectre  are  plotied  over  die  fiequency  band  used  in  the  invereioa  The  spectre  on  the  top  are 
finom  an  earthquake  420  km  fitom  NORESS,  the  spectre  in  the  middle  are  ftom  a  mining  explo> 
Sion  ai  960  km,  and  the  spectre  on  the  bottom  are  from  an  earthquake  at  1 170  km.  The  qiec- 
tral  amplitudes  are  in  nanometer-seconds. 


Table  2.2.  Inversion  Results 


Q(f) 

Geometric  Spreading 

Source  Terms 

Qo 

n 

ro 

m 

log  AfS*P 

K  log  Afg’ 

Pn: 

325 

1.0 

B 

1.35A/i  +  16.9t 

0.26 

100.0 

B 

1.35Ml  +  16.9t 

0.27  l.Q3ML+n.l 

t  Fit  by  eye. 


2.3.2  Simulation  of  NORESS  spectra 

Based  on  the  inversion  parameters  in  Table  2.2,  it  is  possible  to  simulate  spectra  at 
NORESS  for  a  chosen  source  and  epicentral  distance.  For  example.  Figure  2.4  displays 
predicted  Pn  and  Lg  spectra  at  distances  of  SOO,  800.  and  1000  km  for  an  =  3.0  explosion. 
The  lower  curve  is  the  average  ambient  noise  spectrum  at  NORESS  [Suteau-Henson  and 
Bache,  1988].  The  Pn  spectrum  is  nearly  parallel  to  the  noise  curve  at  range  of  about  SOO  km. 
This  is  consistent  with  a  NORSAR  study  that  examined  data  from  the  high-frequency  element 
of  the  NORESS  array  [Ringdai  et  al.,  1986].  At  longer  ranges  the  frequency  of  the  maximum 
signal-to-noise  ratio  (SNR)  is  considerably  less  than  IS  Hz.  For  example,  for  distances  greater 
than  KXX)  km  the  predicted  frequency  of  the  maximum  SNR  is  between  4-5  Hz.  Thus  our 
results  for  NORESS  do  not  support  the  main  conclusion  of  Evernden  et  al.  [1986].  that  Pn 
detection  capability  would  be  enhanced  at  frequencies  greater  than  20  Hz.  except  at  distances 
less  than  400-5(X)  km. 

A  similar  interpretation  for  Lg  is  not  possible  since  we  don’t  have  a  parameterization  of 
pre-Lg  noise  (the  ambient  noise  curve  is  plotted  with  the  Lg  spectra  in  Figure  2.4  only  for 
comparative  purposes).  The  pn-Lg  noise  is  non-stationaiy  because  it  includes  the  coda  of  pre¬ 
viously  arriving  phases.  Here  we  simply  note  that  the  Lg  spectrum  reaches  the  noise  level  at  a 
much  lower  frequency  than  the  Pn  spectrum,  while  at  long  periods  the  Lg  amplitude  exceeds 
that  of  Pn  by  as  much  as  a  factor  of  10.  This  is  consistent  with  the  observation  that  regional 
seismograms  are  characterized  by  Lg  being  the  largest-amplitude  phase  and  Pn  having  a  higher 
domiruuit  frequency. 
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Predicted  Pn  Spectra  Predicted  Lg  Spectra 
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2.4.  Predicted  NORESS  (8)  Pt$  wnil(b)Lg  displacement  spectn  at  three  epicentral  dis- 
for  a  magnitude  3.0  explosion,  based  on  the  inversion  results.  The  average  NORESS 
ivas  estimated  from  43  5-s  samples  taken  prior  to  Pn  [Sutem-Henson  and  Bache,  1988], 


3.  TEMPORAL  AMPLITUDES  AND  DETECTABIUTY 


Signal  detecticm  at  NORESS  is  declared  when  the  short  term  average  (ST A)  divided  by 
the  long  teirn  average  (LTA)  exceeds  a  predetennined  threshold.  These  amplitudes  are  meas¬ 
ured  on  filtered  beams.  Therefore,  to  use  our  spectral  parameterization  to  predict  and  extrapo¬ 
late  NORESS  detection  capability,  it  is  necessary  to  determine  the  relationship  between  single¬ 
channel  spectral  amplitudes  and  the  time-domain  amplitudes  used  in  signal  detection.  In  this 
section  we  determine  this  relationship  and  use  it  to  compare  predicdons  based  on  our  inversion 
results  to  observed  NORESS  detection  capability. 

We  start  by  comparing  observed  detection  statistics  at  NORESS  to  predictions  based  on 
our  spectral  parameterization.  Next,  we  determine  the  relationship  between  temporal  and  spec¬ 
tral  amplitudes  on  filtered  beams.  In  Section  3.2.3,  we  combine  this  with  estimates  of  the 
fiequency-depettdent  i-eam  gain  to  develop  an  expression  for  the  temporal  SNR  on  a  filtered 
beam  in  terms  of  single-channel  spectral  amplitudes.  In  Section  3.3,  we  define  a  function 
called  the  "detectability"  that  gives  simple  and  accurate  estimates  of  the  NORESS  detection 
capability  directly  from  the  observed  temporal  SNR  (=STA/LTA). 

3.1  Detection  statistics 

Detection  and  post-detection  processing  of  NORESS  data  is  done  by  the  automated 
array-processing  package,  SAIAP.  The  detection  processing  is  identical  to  that  in  RONAPP 
[Mykkeltveit  and  Bungum,  1984],  and  uses  a  conventional  STA/LTA  detection  algorithm 
applied  to  a  set  of  filtered  beams.  The  output  of  SAIAP  includes  the  STA  and  LTA  for  each 
detection  measured  on  the  standard  beam  with  greatest  signal-to-noise  ratio  (SNR).  The 
NORESS  array  configuration  is  shown  in  Figure  3.1.  Table  3.1  lists  the  standard  beams  used 
by  SAIAP.  Beams  1-7  and  17-20  are  infinite-velocity  beams,  8-16  are  steered  beams,  and  21- 
24  are  formed  from  horizontal  components.  No  more  than  three  steering  azimuths  are  used  for 
a  given  frequency  filter.  Beams  1-17  are  coherent  beams  and  18-24  are  incoherent  beams. 
Coherent  beams  are  formed  by  delaying  each  channel  by  the  proper  amount  determined  from 
the  steering  azimuth  and  velocity,  summing  all  channels,  and  band-pass  filtering.  Incoherent 
beams  are  fotmed  by  band-passing  each  channel,  delaying  and  rectifying,  and  sununing  all 
channels.  Incoherent  beams  are  particularly  well-suited  for  detection  of  signals  with  low 
coherency  across  the  array  [Ringdal,  198Sb].  For  example,  all  but  one  of  the  Pn  detections 
used  in  this  study  are  on  coherent  beams,  whereas  68%  of  the  Lg  detections  are  on  incoherent 
beams. 

The  most  fundamental  observation  that  our  spectral  parameterization  should  reproduce  is 
the  distance  dependence  of  the  frequency  of  the  maximum  SNR.  This  frequency  depends  pri¬ 
marily  on  the  Q,  since  most  of  the  events  have  ^  3.0  and  high  apparent  comer  frequen¬ 
cies.  Figure  3.2  shows  the  percentage  of  Pn  and  Lg  detections  for  several  frequency  bands  at 
four  ranges.  Although  a  phase  may  be  detected  on  many  beams,  it  is  plotted  in  Figure  3.2 
only  for  the  hequetKy  band  for  the  beam  with  the  maximum  SNR.  As  expected,  detections 
from  close  events  generally  occur  on  a  higher-frequency  beam  than  those  from  events  at  larger 
distances.  For  example,  the  maximum  SNR  for  Pn  typically  occurs  at  a  frequency  greater  than 
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Figure  3.1.  N(MtESS  amy  conflguntioii.  The  irrey  coosisa  of  25  shoit-period  umrumena  in 
concentric  ring!  widi  a  maximum  diameter  of  3  km.  The  array  was  designed  for  the  edunoe- 
ment  of  fiequencies  between  U  and  5.0  Hz  [MyUoeltveit,  1983]. 
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TaUe  3.1.  Standaid  Beams  used  by  SAIAP.  Listed  are  beam  number,  steering  azimuth  (deg) 
and  slowness  (Vkm),  frequency  filter  (Hz),  beam  type  (coherent  or  incoherent),  sensor  weights, 
and  detector  thresholds.  The  individual  sensors  are  identified  by  ring  (A-D),  sensor  number 
(0-9),  and  component  (Z.  N.  or  E).  The  thresholds  were  determined  from  false  alann  statistics 
[Kvaema  et  al.,  1987a]. 
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of  with  the  maximum  SNR  for  each  ftequency  filter  grouped 

e$  for  (a)  Pti  and  (b)  Lg.  The  disumce  imeivals  and  the  number  of 
are  indicated  on  the  right 


8  Hz  at  distances  less  than  400  km  and  between  3-S  Hz  in  the  700-1000  km  distance  range. 
This  is  consistent  with  the  results  of  Ringdal  [198Sa]  who  found  that  the  best  SNR  for  Pn 
increased  from  3-5  Hz  at  about  1000  km  to  more  than  8  Hz  at  local  distances.  Note  that  the 
predicted  Pn  spectra  shown  in  Figure  2.4  are  generally  consistent  with  these  observations.  A 
similar  comparison  for  Lg  is  not  possible  since  we  don’t  have  a  parameterization  of  the  pre-Lg 
noise  spectrum.  However  there  is  at  least  qualitative  agreement  between  the  predicted  and 
observed  dominant  frequency  of  Lg'. 

32  Relationship  between  spectral  and  temporal  amplitudes 

In  this  section  we  determine  the  relationship  between  temporal  and  spectral  amplitudes  on 
filtered  beams,  and  combine  this  with  estimates  of  the  frequency-dependent  beam  gain  to  deter¬ 
mine  die  relationship  between  the  SNR  used  in  signal  detection  and  the  single  channel  spectral 
SNR.  Since  Lg  is  strongly  dispersed  and  the  pre-Lg  LTA  is  non-stationary,  we  are  unable  to 
parameterize  the  Lg  SNR  in  a  form  suitable  for  extrapolation.  Therefore,  in  this  section  we 
concentrate  on  Pn. 

32.1  Noise 

Qualitativdy,  the  LTA  is  a  measure  of  the  RMS  amplitude  on  a  filtered  beam  averaged 
over  about  30  s  prior  to  a  detection.  The  LTA  is  updated  every  0.S  s  and  is  expressed  as  a 
weighted  sum  of  past  ST  As.  There  is  a  simple  relationship  between  the  RMS  amplitude  of 
random  noise  and  its  power  spectral  density  [e.g.  Aki  and  Richards,  1980;  p.  493].  From 
Paiseval’s  theorem 


Akms  *  (2  PSD  Af)*'*  (3.1) 

where  is  the  bandwidth  and  PSD  is  the  power  spectral  density.  This  equation  assumes  that 
the  PSD  is  fiat  over  the  bandwidth  Af.  The  LTA  for  Pn  is  an  approximate  measure  of  the 
RMS  ambient  noise  amplitude  on  a  filtered  beam  and  should  be  approximately  related  to  the 
PSD  of  the  beam  by  (3.1),  wdiere  Af  is  the  filler  bandwidfii. 

To  test  the  :q)plicability  of  (3.1)  to  pre-Pn  (or  ambient)  noise,  we  computed  ^ctra  of 
unfiltered  beams  over  various  sub-arrays.  Figure  3.3  compares  the  average  noise  spectra 
recorded  on  unfiltered  beams  over  the  four  sub-arrays  used  in  the  standard  beams  to  the  aver¬ 
age  single-channel  noise  spectrum  [Suteau-Henson  and  Bache,  1988].  Fbr  frequoicies  less 
than  about  4  Hz  the  noise  level  is  significantly  reduced  by  including  the  two  outer  rings  of  the 
array.  In  particular,  the  beam  including  sensors  on  only  the  C  and  D  rings  gives  better  noise 
suppression  than  can  be  achieved  using  the  whole  array  [Kvaerna  and  Mykkeltveit.  1986]. 
However,  beyond  5-6  Hz  only  marginal  improvemem  in  noise  suppression  can  be  achieved  by 
including  the  outer  rings.  Superimposed  on  the  spectra  are  the  squared  LTAs  divided  by  2^ 
for  the  186  Pn  detections.  The  LTAs  are  corrected  for  the  instrument  response  and  plotted  at 
the  frequency  of  the  mean  power  in  the  bandwidth  df.  Note  that  the  LTAs  measured  on  the 
8-16  Hz  beam  are  consistent  with  (3.1)  at  a  frequency  of  about  9.5  Hz.  This  is  an  important 
distinction  for  detection  capability  assessment  This  implies  that  the  temporal  SNR  cannot  be 


25 


10® 

10® 

10^ 


Noise  Spectra  and  LTA 


Figure  3.3.  Avenge  NORESS  ambiem  noise  power  spectnl  density.  The  array-avenged 
singk-channel  spectrum  is  fiom  Suteau-Henson  and  Bache  [1988].  The  beam  spectn  are  aver¬ 
ages  of  10  nc^  samples  calculated  for  the  various  sub-amys.  AU  spectn  are  comcted  for 
the  insmiment  response.  The  symbols  are  Pit-LTAs  divided  by  2Af  for  186  events. 


predicted  based  on  a  spectral  parameterization  of  the  signal  and  the  noise  at  the  same  fre¬ 
quency,  unless  narrow  band  filteis  are  used  in  beamforming.  Otherwise  the  dominant  signal 
frequency  may  exceed  the  dominant  noise  frequency  within  the  filter  bandwidth.  This  is  dis¬ 
cussed  in  more  detail  in  a  later  section. 

The  noise  suppression  is  defined  as  the  ratio  of  the  beam  power  spectrum  to  the  array- 
averaged  single-channel  power  spectrum  [Fyen,  1986].  In  general,  we  find  that  at  least  1/N 
noise  suppression  (where  Af  is  the  number  of  array  elements  used  in  beamfotming)  is  achieved 
over  the  frequency  band  used  for  each  of  the  sub-arrays.  However,  the  signals  are  also 
suppressed  by  beamfotmiitg,  and  this  must  also  be  considered  in  establishing  the  frequency- 
dependent  beam  gain  (Section  3.2.3). 

The  pre-Lg  noise  consists  of  ambient  noise,  P  coda,  and  Sn  coda.  Therefore,  we  do  not 
expect  a  simile  relationship  between  Lg-LTA  aral  the  ambient  noise  spectrum.  Furthermore, 
we  don’t  have  a  parameterization  of  the  single-channel  pre-Lg  itoise  spectrum,  so  there  is  no 
reason  to  relate  temporal  and  spectral  amplitudes  as  was  done  for  Pn.  Therefore,  we  simply 
parameterize  the  Lg-LTA  in  terms  of  distance,  frequency,  magnitude,  and  beam  type  [Sereno 
and  Bratt,  1988].  We  find  that  the  Lg-LTA  decreases  with  increasing  distance  and  frequency, 
but  that  the  most  important  effect  for  detection  capability  assessment  is  that  of  magnitude. 
Figure  3.4  [dots  Lg-LTA  versus  NORESS  kx;ai  magnitude.  The  best-fitting  straight  line  to  log 
LTA  as  a  function  of  magnitude  has  a  slope  near  0.8  for  both  coherent  and  incoherent  beams. 
Since  both  Lg  signal  and  "noise*'  increase  widi  increasing  M^,  this  implies  that  the  detectability 
of  Lg  is  much  less  dependent  on  magnitude  dian  that  of  Pn. 
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322  Signal 

The  Pn  spectra  used  in  this  study  ate  those  automatically  computed  by  SAJAP  for  each 
detected  phase.  They  are  computed  for  a  S>s  window  starting  0.3  s  before  the  detection  time 
on  the  vertical  componoit  and  are  noise-corrected  and  averaged  across  the  array.  In  this  sec¬ 
tion,  we  examine  the  relationship  between  these  spectral  amplitudes  and  the  time-domain  STA 
amplitude  on  a  filtered  beam.  However,  the  Lg  energy  distributed  over  the  fixed  group  velo¬ 
city  wirxlow  used  in  the  inversion  is  not  an  appropriate  measure  for  detection  capability  assess¬ 
ment  We  find  that  the  Lg  spectral  amplitudes  computed  for  S-s  windows  are  more  appropriate 
for  that  api^ication.  Specifically,  we  find  that  the  Lg-STA  are  proportional  to  these  spectral 
amplitudes,  and  that  the  proportionality  constant  depends  only  on  beam  type  [Serena  and  Bratt, 
1988].  However,  since  Lg  Q  estimates  based  on  these  spectra  are  not  available  for  other 
regions,  our  results  for  Lg  are  not  eaaly  extrapolated.  Therefore,  the  remainder  of  this  section 
concentrates  on  the  relationship  between  spectral  and  temporal  Pn  amplitudes. 

The  STA  is  defined  as  the  average  absolute  amplitude  in  a  1-s  window  measured  on  the 
beam  with  the  maximum  SNR  and  is  updated  every  sample.  As  a  tough  approximation,  the 
amplitude  of  a  wavelet  is  the  product  of  the  tunplitude  spectral  density  and  the  bandwidfit  [Aid 
and  Richards,  1980;  p.  492].  Therefore,  we  express  the  STA  in  terms  of  the  amplitude  spec¬ 
tral  density  of  the  beam  as 


57-A=A^o)Ar5(/b) 


(3.2) 


where  /o  is  the  dominant  signal  frequency,  is  the  filter  baitdwidth,  A^  is  the  amplitude  ^c- 
trum  of  the  beam,  and  5  is  a  correction  term  to  accoum  for  dispersion  and/or  scattering.  If  all 
of  the  energy  in  the  bandwidth  arrives  at  foe  same  time  then  5-1,  otherwise  5  <  1.  If  the 
signal  is  purely  random,  then  5  can  be  found  from  (3.1).  The  range 


^  5  ^  1.0 


(3.3) 


where  T  is  the  noise  window  length  On  this  case,  5  s),  corresponds  to  that  between  a  purely 
random  and  perfectly  coherent  Pn  signal.  The  lower  bound  on  5  is  0.4S,  0.31,  0.22,  for  2,  4, 
and  8  Hz  bandwidths.  respectively.  Figure  35  plots  examples  of  the  STA  corrected  for  the 
bandwidth  and  5.  superimposed  on  the  spectra  at  the  frequency  of  the  mean  signal  errergy  in 
foe  bandwidth  AT-  Comparing  the  STAs  to  the  beam  spectra  for  10  events,  we  find  that  a  con¬ 
sistent  value  of  5  =  0.7  relates  Pn  temporal  and  spectral  amplitudes  on  beams  filtered  between 
2  and  5  Hz.  However,  between  8-16  Hz  we  find  5  ==  0.3.  This  lower  value  is  probably 
caused  by  the  rapid  decay  of  the  beam  spectra  in  this  frequency  band,  which  reduces  the 
"effective"  signal  bandwidth.  Increased  scaoeiing  at  high  frequency  would  also  have  the  effea 
of  lowering  5. 
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Figure  3.5.  Pn  beam  spectra  for  10  events  recorded  at  NORESS.  The  five  spectra  on  top  are 
for  beams  filtered  between  approximately  2  and  5  Hz.  The  lower  spectra  are  for  hMtn«  filtered 
between  8  and  16  Hz.  The  asterisk  on  each  ^ctnim  |dots  STA/(AA5)  where  5  «  0.7  between 
2-5  Hz  and  5  «  0.3  for  the  8-16  Hz  beam. 


32  J  Signal-to-Noise 

A  major  objective  of  this  study  is  to  predict  the  SNR  as  a  function  of  distance,  magnitude 
and  frequency  based  on  a  parameterization  of  array-averaged  spectra.  To  do  this,  we  must 
address  two  separate  issues.  One  is  the  relation  between  temporal  and  spectral  amplitudes  on 
filtered  beams,  and  the  other  is  the  relation  between  beam  spectra  and  single-channel  spectra. 
In  Sections  3.2.1  and  3.2.2  we  addressed  the  first  issue  by  expressing  the  Pn  LTAs  and  STAs 
in  terms  of  the  beam  spectral  amplitudes.  In  this  section  we  investigate  the  second  issue  which 
involves  the  determination  of  the  frequency-dependent  beam  gaia 

While  beamfoiming  can  produce  l/N  or  greater  noise  suppression,  it  also  suppresses 
uncoirelated  high-frequency  signal  energy.  The  beam  gain  is  defined  as  the  ratio  of  the  SNR 
on  the  beam  and  the  single-channel  SNR.  Therefore,  using  (3.1)  and  (3.2)  we  can  express  the 
SNR  on  a  filtered  beam  in  terms  of  the  array-averaged  single-channel  spectra  as 

where  A(f)  and  N(f)  are  single-channel  signal  and  noise  spectra  and  G(f)  is  the  frequency- 
dependent  beam  gairt  The  only  term  in  (3.4)  that  we  have  not  estimated  is  the  beam  gaia 
Frinn  the  observed  SNR  and  our  estimates  of  5.  we  can  determitte  the  average  G(f)  for  specific 
beams  using  (3.4).  Note  that  the  signal  and  rtoise  are  not,  in  gerteral,  evaluated  at  the  same 
frrequency  (each  is  evaluated  at  the  frequency  of  the  mean  ertergy  or  power  over  the  bandwidth 
AO-  Figure  3.6  plots  CUf)  for  the  four  combinations  of  sub-array  and  filter  bandwiddt  typical 
for  Pn  detections.  The  mean  gain  for  each  beam  type  is  shown  as  die  solid  horizontal  line. 
Table  3.2  summarizes  the  beam  gain  results. 


Table  3.2.  Beam  gairt 


Sub-array 

Beam  Numbers 

FrequetKy 

AT 

Detections 

<N 

<G> 

o 

5 

AO-BCD 

3,8-16 

2.0 

68 

4.7 

3.0 

1.4 

0.7 

AO-BC 

4-5 

2.0 

27 

3.6 

3.2 

2.6 

0.7 

AO-AB 

6 

4.0-8.0 

4.0 

30 

3.0 

1.11 

0.8 

0.7 

AO-AB 

7 

8.0-16.0 

8.0 

56 

3.0 

1.09 

0.6 

0.3 

For  beams  with  center  frequencies  below  5  Hz,  the  beam  gain  approadies  although 
there  is  much  scatter.  The  noise  suppression  is  much  more  consistent  about  V^.  The  large 
scatter  in  Figure  3.6  is  due  to  variability  in  signal  coherence.  However,  the  gains  at  higher  fre¬ 
quency  are  near  unity.  That  is,  the  array  offers  little  or  no  advantage  over  single  channels. 
This  is  consistent  with  regional  beam  gain  spectra  displayed  by  Mykkeltveit  et  al.  [1985]. 
Since  the  gains  listed  in  Table  3.2  depend  on  our  value  of  5,  this  consistency  provides 
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Figure  3.6.  Beam  gain  calculated  using  (3.4)  for  the  four  combinadons  of  sub-amy  and  fre¬ 
quency  filter  typical  for  Pn  detection.  The  means  are  plotted  bounded  by  one  standard  devia¬ 
tion.  A  bold  line  at  is  ironed  for  reference.  Note  that  the  same  axes  scaling  is  used  for  all 
{dots  and  that  not  all  of  die  points  used  to  detennine  the  mean  and  standard  deviation  are 


independent  suppoit  for  our  estimate  of  that  parameter.  The  NORESS  array  was  designed  to 
enhance  signal  frequencies  in  the  range  1.5-5.0  Hz  [Mykkeltveit,  1983].  Thus,  it  is  not  surpris¬ 
ing  that  the  array  does  not  offer  enhanced  SNR  at  higher  frequencies.  This  simply  means  that 
the  signal  and  the  noise  are  both  uncorrelated  over  the  sub-array  used  in  beamforming.  Note, 
however,  that  it  is  likely  that  the  high-frequency  beam  gain  could  be  improved  by  adding 
steered  beams  and/or  by  filtering  over  a  narrower  frequency  band.  This  is  discussed  in  Section 
3.4. 

3  J  Regional  wave  detectability 

In  this  section  we  introduce  an  empirical  parameterization  of  the  temporal  SNR  that  we 
call  the  "detectability,"  which  gives  the  probability  of  detection  of  a  given  magnitude  event  as 
a  function  of  epicentral  distance.  The  detectability  is  computed  directly  from  parameters  that 
defitK  the  way  an  automatic  detector  works  (e.g..  the  ST A,  LTA,  and  pre-set  detector  thres¬ 
holds).  Therefore,  it  is  possible  to  very  accurately  simulate  the  detection  capability  of  arrays 
equipped  with  automatic  array-processing  software.  However,  another  question  is  how  to 
extrapolate  these  results  to  areas  with  different  attenuation.  For  this  we  use  the  relationship 
between  the  temporal  SNR  and  the  spectral  inversion  parameters  derived  in  the  previous  sec¬ 
tion.  To  validate  this  procedure,  we  compare  "predicted”  (based  on  our  inversion  results)  and 
"observed"  Pn  detectability  at  NORESS.  Since  we  don’t  have  a  similar  relationship  for  Lg,  we 
simply  present  an  empirical  parameterization  of  NORESS  Lg  detectability  and  note  that  the 
results  cannot  be  extrapolated  with  confidence  to  other  regions. 

The  automatic  detection  of  a  seismic  signal  occurs  whcm  the  STA/LTA  exceeds  a 
predetermined  threshold.  Therefore,  dividing  by  the  threshold  provides  a  cmnmon  basis  for 
comparing  the  SNR  measured  on  different  beams  [Kvaerna  et  ai,  1987a].  Since  the  beam 
thresholds  are  determined  from  false  alarm  statistics,  this  means  that  the  probability  of  a  spuri¬ 
ous  detection  is  the  same  on  all  beams.  We  define  the  detectability,  D(A),  as 


D(A)  =  log 


SNR 

th 


-olMl 


(3.5) 


where  SNR  is  the  maximum  STA/LTA  on  all  detecting  beams,  th  is  the  beam  threshold,  and  a 
is  chosen  such  that  Z)(A)  does  not  depend  on  the  source.  The  value  of  a  is  close  to  one  for  Pn 
but  much  less  than  one  for  Lg  because  both  the  STA  and  the  LTA  depend  on  Ml..  An  impor¬ 
tant  aspect  of  the  detectability  is  that  it  includes  frequency  dependence.  Fbr  example,  events 
within  400-5(X)  km  of  NORESS  typically  have  the  maximum  SNR  for  Pn  on  a  high-frequency 
beam,  while  more  distant  events  are  usually  detected  with  the  maximum  SNR  on  a  lower  fre¬ 
quency  beam  (Figure  3.2).  The  detectability  can  be  interpreted  simply  in  terms  of  the  magni¬ 
tude  detection  threshold.  For  example,  when  the  SNR  is  equal  to  the  beam  threshold,  -£>(A)  is 
equal  to  a  multiplied  by  the  50%  NORESS  Ml  detection  threshold.  Assuming  that  both  signal 
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and  noise  are  log  normally  distributed,  the  probability  of  detecting  wave  k  from  source  j  is 


PkiTj)  O-k  Mlj 

o* 


(3.6) 


where  (t>(x)  is  the  Gaussian  cumulative  distribution  function  [e.g.,  Abramowitz  and  Scgun, 
1964;  p.  931]  and  o*  is  the  standard  deviation  of  the  detectability  function.  This  relationship  is 
a  valid  representation  of  the  probability  of  detection  below  the  source  comer  frequency. 

33.1  Pn  Detectability 

Figure  3.7  plots  the  temporal  Pn  detectability  defined  by  (3.5)  with  a  =  1.0.  Only  explo¬ 
sions  are  used  because  the  earthquake  detectability  can  be  biased  by  the  unknown  radiation  pat¬ 
tern.  That  is,  it  is  not  possible  to  deteimiire  an  a  for  earthquakes  that  effectively  removes  the 
source  contribution  to  tire  SNR.  We  also  excluded  events  with  ^  3.0  that  were  detected  on 
one  of  the  two  highest-frequency  beams  because  their  corner  frequencies  are  expected  to  be 
within  the  filter  bandwidth.  The  logarithmic  decay  of  the  Pn  detectability  is  approximated  by 
D(A)  =  -0.S3  -  1.93  logA,  where  A  is  epicentral  distance  in  degrees  and  the  coefficients  were 
obtained  by  least-squares.  This  is  plotted  in  Figure  3.7  bounded  by  one  standard  deviation  (o 
=  0.31).  Based  on  this  parameterization,  tiie  50%  NORESS  Mi,  threshold  for  detecting  Pn  at 
400  km  is  about  1.6  and  at  1000  km  is  about  2.4. 

In  order  to  validate  the  relationship  between  the  temporal  and  spectral  SNR  derived  in  the 
previous  section,  we  compare  the  [xedicted  (based  on  the  spectral  inversion  results)  and 
observed  NORESS  Pn  detectability.  We  use  our  parameterization  of  the  single-channel  spec¬ 
trum  A(f),  and  the  average  NORESS  noise  spectnun  N(f)  (Figure  2.4),  to  predict  the  detection 
capability  of  specific  beams  using  (3.4).  We  use  gain  for  beams  with  frequency  filters 
between  2-5  Hz  and  a  gain  of  1.0  for  the  two  highest  hequency  beams.  We  set  a  =  1.35, 
which  is  the  slope  of  the  log  M(fMi  relation  for  explosions  (Figure  2.2).  Note  that  the 
difference  between  the  implied  50%  NORESS  Mi,  threshold  for  a  =  1.35  and  a  =  1.0  is  only 
about  0.1.  Figure  3.8  demonstrates  the  general  consistency  between  the  predicted  and  observed 
Pn  detectability. 

Although  the  hequency  of  the  maximum  SNR  predicted  from  the  spectral  inversion  is 
greater  than  10  Hz  for  distances  less  than  800  km  (Figure  2.4a),  the  predicted  frequency  of 
optimal  detectability  is  much  less  than  this.  This  is  because  the  single-channel  SNR  advantage 
at  high  frequency  is  erased  by  a  lower  beam  gain.  As  a  result,  our  model  predicts  greater 
detection  capability  on  2.5-5  Hz  beams  than  on  the  higher-frequency  beams  for  ranges  greater 
than  about  200  km.  However,  the  observatitms  show  a  significant  number  of  Pn  detections  on 
the  higher-frequency  beams  at  longer  ranges,  even  though  their  SNRs  are  consistent  with  the 
model.  A  likely  explanation  is  that  spectral  modulation  caused  by  ripple-firing  with  delays 
between  150-250  ms  consistently  suppresses  amplitudes  less  than  4-5  Hz  for  mine  blasts 
[Baumgardt  and  Zeigler,  1987].  In  many  cases  our  Q  model  correctly  predicts  the  hi^  fre¬ 
quency  (>  5  Hz)  Pn  amplitude  of  these  events,  but  overestimates  the  spectral  level  between  2-5 
Hz  [Sereno  and  Bratt,  1988].  Since  the  time  lags  that  produce  4-5  Hz  spectral  modulations  are 
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Figure  3.7.  Pn  detectability  defined  by  (3.S)  with  a  »  1.0  for  102  exjdosions.  Seven  events 
were  excluded  with  ^  3.0  because  their  comer  frequencies  are  expected  to  be  within  the 
filter  bandwidth.  The  curve  is  the  logarithmic  parameterizatioi)  of  D(A)  bounded  by  oik  stan¬ 
dard  deviation. 
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Cpicentral  Olstane*  (degrtas) 


Hgufe  3.8.  Comparison  of  predicted  (based  on  our  spectnd  inversion)  and  observed  Pn  detec¬ 
tability  defined  by  (3.5)  with  a  «  1.35.  This  value  of  a  is  based  on  the  observed  rate  of 
increase  of  exfdoskxi  momeru  with  M/,  (see  Hgure  2.2). 
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consistent  with  the  destnictive  interference  of  multiple  explosions,  we  interpret  this  as  an 
unmodeled  source  complexity  rather  than  a  path  effect. 

3J.2  Lg  Detectability 

In  this  section  we  present  the  results  of  an  empirical  parameterization  of  NORESS  Lg 
detectability.  However,  without  a  parameterization  of  the  pre-Lg  noise  spectrum  Gike  Figure 

3.3  for  Pn),  we  caiuiot  predict  the  Lg  detectabOity  based  on  Q  and  moment.  Therefore,  while 
our  results  accurately  represent  Lg  detectability  at  NORESS,  they  cannot  be  extrapolated  with 
confidence  to  other  regions.  The  value  of  a  in  (3.S)  that  effectively  removes  the  source  contri¬ 
bution  to  the  detectability  is  equal  to  the  difference  in  slopes  of  the  log  moment-magnitude 
(~  1.0)  and  log  LTA-magnitude  relations  C0.8).  Therefore  we  expect  the  appropriate  value  of 
a  to  be  close  to  0.2.  If  the  source  contribution  is  effectively  removed,  the  detectability  should 
not  depend  on  magnitude,  provided  the  measurements  are  taken  below  the  comer  frequency. 
We  find  that  subtracting  values  much  greater  than  0.2  produces  an  obvious  magiutude 
dependence  in  the  Lg  detectability,  indicating  that  the  source  contribution  to  the  SNR  was  not 
adequately  removed  [Sereno  and  Bratt,  1988].  Figure  3.9  plots  the  observed  Lg  detectability 
defined  by  (3.S)  with  a  =  0.2.  We  include  both  earthquakes  and  explosions  since  Lg  is  rela¬ 
tively  insensitive  to  the  source  radiation  pattern.  The  logarithmic  decay  of  the  Lg  detectability 
with  a  =  0.2  is  approximated  by  D(A)  =  -0.11  -  0.64  log  A,  where  the  coefficients  were 
obtained  by  least-squares.  This  is  plotted  in  Figure  3.9  bounded  by  one  standard  deviation  (o 
=  0. 12).  From  this  parameterization,  the  50%  NORESS  A//,  threshold  for  detecting  Lg  is  about 

1.2  at  400  km  and  about  2.S  at  1000  km. 

3.4  Detection  capability  of  the  NORESS  array 

In  this  section  we  estimate  the  detection  capability  of  NORESS  using  our  parameteriza- 
tioa  These  estimates  are  compared  to  the  results  of  a  NORSAR  study  that  estimated  detection 
capability  by  comparing  the  NORESS  detections  to  bulletins  produced  by  local  seismic  net¬ 
works  in  Feruioscandia. 

It  is  straightforward  to  use  our  parameterization  of  0(A)  to  estimate  the  detection  capabil¬ 
ity  of  the  NORESS  array.  We  use  (3.6)  together  with  the  detectability  to  calculate  die  proba¬ 
bility  of  detecting  a  given  magnitude  event  as  a  function  of  epicentral  distance.  For  example, 
the  90%  NORESS  detection  threshold  for  Pn  based  on  D(A)  with  a  =  1.0  is  2.0  at  4(X)  km 
and  2.8  at  1000  km.  For  Lg  the  90%  detection  threshold  with  a  =  0.2  is  2.0  at  400  km  and 

3.3  at  1000  km.  The  90%  thresholds  for  detection  of  either  Pn  or  Lg  are  1.7  at  4(X)  km  and 
2.6  at  10(X)  km.  These  estimates  are  generally  consistent  with  those  of  Ringdal  [1986]  who 
estimated  the  90%  threshold  for  detecting  P  phases  between  700  and  14(X)  km  (average 
about  960  km)  to  be  2.7  and  for  detecting  P  or  secondary  phases  in  the  same  distance  range  to 
be  2.5.  We  have  not  studied  the  detectability  of  Sn,  so  we  carmot  directly  compare  our  secon¬ 
dary  phase  results  to  his.  Since  the  addition  of  Sn  can  only  improve  the  detectability,  our 
results  for  Pn  or  Lg  detection  give  an  upper  bound.  Table  3.3  compares  the  results  of  the  two 
studies. 
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Table  3.3.  Regional  wave  detection  capability  at  NORESS. 


Phase 

50%  Threshold 

this  study  Rtngdo/ [1986] 

90%  Threshold 

this  study  Ringdal  [1986] 

P  «ily 

2.4 

2.3 

2.8 

2.7 

Lg  only 

P  or  Secondary 

2.5 

3.3 

Phases 

2.2t 

1.9 

2.6t 

2.5 

t  This  study  uses  Lg  as  the  only  secondary  phase,  Ringdal  [1986]  used  Sn  or  Lg. 


In  Figure  3.10  the  detection  capability  estimates  based  on  our  temporal  parameterization 
are  compared  to  those  based  on  the  spectral  inversion.  The  dashed  curves  are  our  estimates  of 
the  90%  NORESS  thresholds  for  detecting  Pn  and  for  detecting  Pn  or  Lg,  based  on  the 
temporal  detectability  curves  in  Figures  3.7  and  3.9.  The  solid  curve  is  the  90%  threshold 
for  detecting  Pn  based  on  the  detectability  predicted  by  the  spectral  inversion  results  (Figure 
3.8).  Also  indicated  in  Figure  3.10  are  the  results  of  Ringdal  [1986].  In  general,  these  three 
independent  studies  produce  consistent  estimates  of  the  regional  wave  detection  capalnlity  of 
the  NORESS  array. 

Note  that  our  estimates  of  the  detection  capability  of  the  NORESS  array  are  based  on  the 
standard  beam  set  listed  in  Table  3.1.  However,  it  is  possible  that  NORESS  capability  could 
be  enhanced  by  including  more  steered  beams  [T.  Kvaema,  personal  communication,  1988],  or 
by  modifying  the  frequency  filters  [Kvaerna  et  al.,  1987b;  Sereno  and  Bran,  1988].  For  exam¬ 
ple,  the  predicted  Pn  spectrum  at  ranges  less  than  SOO  km  decays  less  rapidly  than  the  noise 
spectrum  beyond  8  Hz  (Figure  2.4).  This  suggests  that  the  Pn  detection  capability  of  the  8-16 
Hz  beam  (beam  7,  Table  3.1)  could  be  improved  by  filtering  at  higher  frequency.  This  is 
confirmed  by  Kvaerna  et  al.  [1987b]  who  compiled  detection  statistics  for  events  in  western 
Norway.  Most  of  their  events  were  at  ranges  less  than  450  km.  They  found  that  60%  of  their 
Pn  detections  occurred  with  the  maximum  SNR  on  a  10-16  Hz  beam  while  only  4%  occurred 
on  the  8-16  Hz  beam.  Therefore,  our  results  may  underestimate  the  detection  capability  of  a 
NORESS-type  array  with  an  improved  beam  deployment 
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Figure  3.10.  F^mares  of  90%  A#£,  detection  thresholds  at  NORESS  as  a  function  of  epicentral 
distance.  The  dashed  curves  are  based  on  a  parameterization  of  observed  temporal  amplioides. 
Curves  for  detecting  Pn  and  for  detecting  either  Pn  or  Lg  are  (dotted.  The  solid  curve  is  based 
on  the  results  of  our  inversion  of  Pn  spectra.  The  horizontal  bars  are  detection  threshcdds  at 
NORESS  as  detennined  by  com(Mring  NORESS  detections  to  bulletins  (noduced  by  local 
seismic  netwoiks  [Ringdal,  1986].  Events  for  this  study  were  at  ranges  between  700  and  1400 
and  the  average  distance  was  about  960  km  (X). 
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4.  DETECTION  CAPABIUTY  SIMULATION 

In  this  section  we  present  simulations  of  the  detection  capability  of  hypothetical  netwoilcs 
in  the  Soviet  Union.  These  simulations  are  nonnalized  by  the  NORESS  results  {xesented  in 
previous  sections.  Our  parameterization  of  NORESS  detection  capability  in  terms  of  the  noise, 
attenuation,  and  signal  processing  characteristics  allows  us  to  predict  detection  capability  of  a 
NORESS-type  array  for  other  regions  with  different  conditions.  Of  course,  the  validity  of 
these  predictions  depends  on  the  accuracy  of  the  attenuation  and  noise  estimates  for  areas  of 
interest.  Since  these  are  not  available  for  the  Soviet  Union,  we  assume  a  range  of  conditions 
and  determine  the  sensitivity  of  our  results  to  these  assumptions. 

The  network  detection  capability  simulations  presented  in  this  section  were  computed 
using  an  extended  version  of  Seismic  Network  Assessment  Program  for  Detection  [SNAPID, 
Ciervo  et  al.,  198S]  called  SNAPIDX  [Bratt  et  al.,  1987].  The  program  first  calculates  the  pro¬ 
bability  that  a  seismic  i^iase  will  be  detected  at  a  single  station  using  (3.6).  Next,  the  network 
detection  thresholds  are  determined  from  the  detection  probabilities  of  individual  stations.  The 
threshold  for  each  epicenter  in  the  grid  is  determined  by  varying  the  magninide  until  the  detec¬ 
tion  probability  equals  the  desired  confidence  level. 

4.1  Extrapolation  of  NORESS  capability 

Simulations  of  the  detection  capability  of  a  hypothetical  network  of  seismic  stations  in  the 
Soviet  Union,  require  assumptions  about  the  signal  and  noise  characteristics  in  that  tegioa 
The  simplest  assumption  is  to  use  our  results  from  NORESS  for  all  stations.  That  is,  we 
assume  that  all  aspects  of  regional  wave  propagation  (signal  and  noise)  at  all  stations  are  ident¬ 
ical  to  what  we  observe  at  NORESS.  Even  here  we  require  some  extrapolation  since  our 
NORESS  results  are  strictly  valid  only  to  rdiout  12*’.  To  estimate  detectability  at  longer  ranges, 
we  compared  the  Preliminary  Determination  of  Epicenters  (PDE)  bulletin  to  detections  at 
NORESS.  Based  on  observed  arrival  time,  azimuth  and  phase  velocity;  P-wave  detections  at 
NORESS  were  found  for  nine  events  listed  in  the  PDE  bulletin  at  distances  greater  than  12°. 
Table  4.1  lists  the  PDE  locations  and  magnitudes  of  these  events.  Also  listed  are  the  predom- 
irtant  signal  frequencies  determined  by  SAIAP.  Figure  4.1  displays  the  D(A)  calculated  for 
these  events  using  (3.5)  with  a  =  1.0.  These  events  are  considerably  larger  than  those  used  in 
our  regional  study,  and  it  is  possible  that  the  source  corner  frequency  could  bias  the  detectabil¬ 
ity  estimate.  Therefore,  we  assume  corner  frequency  scales  inversely  with  the  cube  toot  of  the 
moment,  and  that  it  equals  10  Hz  for  a  magnitude  3.0  event  (based  on  our  inversion  results),  to 
get  a  rough  estimate  of  the  source  comer  frequencies.  An  oo-square  source  model  was  used  to 
correa  the  amplitudes  for  the  predicted  source  comer  frequency  effect  The  results  are  indi¬ 
cated  by  vertical  bars  on  D(A)  in  Figure  4.1.  The  source-corrected  D(A)  were  parameterized  by 
a  logarithmic  decay  and  the  coefficients  were  determined  by  least-squares.  The  resulting 
teleseismic  P-wave  detectability,  D(A)  =  -0.69  -  1.82  log  A,  is  plotted  as  a  solid  line  in  Figure 
4.1.  Shown  for  comparison  is  the  1-Hz  amplitude-distance  curve  derived  by  Veith  and  Claw¬ 
son  [1972],  normalized  by  the  average  1-Hz  NORESS  noise  (3.1  nm-s),  Vn  beam  gain,  and  the 
filter  bandwidth  and  detector  threshold  used  for  the  NORESS  1-3  Hz  beam.  Note  that  since 
the  predominant  frequency  generally  decreases  with  increasing  epicentral  distance,  the  detecta¬ 
bility  carmot  be  simulated  with  the  amplitude-distance  curve  at  a  single  frequency.  Our  curve 
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Table  4.1.  Eveitts  used  to  estimate  teleseismic  detectabflity 


Event 

Date 

Time 

Location 

Depth 

Range 

m* 

Frequency 

1 

7-18-85 

21:15 

65.97N  40.86E 

0.0 

14.0 

5.0 

2.2 

2 

11-23-85 

20:14 

42.31N  19.88E 

9.0 

19.1 

3.4 

2.5 

3 

8-  1-86 

13:57 

73.04N  5537E 

33.0 

20.5 

4.6 

2.5 

4 

7-25-85 

3:11 

49.89N  78.15E 

0.0 

37.6 

5.0 

2.2 

5 

2-17-86 

12:36 

36.66N  71.17E 

33.0 

44.0 

5.0 

3.1 

6 

10-27-86 

14:12 

7.49N  36.60W 

10.0 

64.1 

5.2 

1.1 

7 

11-15-85 

11:56 

44.15N  148.29E 

33.0 

69.4 

5.2 

3.5 

8 

10-27-86 

12:16 

42.74N  145J5E 

33.0 

70.0 

4.9 

3.6 

9 

8-13-86 

15:27 

7.69N  74.71W 

60.0 

81.5 

5.0 

IJ 
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Talesalsmic  P-Wav«  Oetactabllliy 


Figure  4.1.  Teleseismic  F*wave  detectability.  The  asterisks  are  the  D(A)  computed  using  (3.S) 
with  a  =  1.0  for  the  nine  events  listed  in  Tsfole  4.1.  The  solid  bars  indicate  the  shift  in  D(A) 
after  conecting  for  the  expected  source  comer  frequency.  The  solid  curve  is  the  logarithmic 
paraneteiization  of  the  source-corrected  D(A).  The  triangles  indicate  our  parameterization  of 
the  NORESS  regional  Pn  detectability  shown  in  Figure  3.7.  The  squares  indicate  1-Hz  detec¬ 
tability.  based  on  the  amplitude-distance  curve  of  Veith  and  Clawson  [1972],  NORESS  1-Hz 
noise,  a  2  Hz  filter  bandwidth,  and  Vn  beam  gaia 
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indicates  larger  amplitudes  (or  greater  detection  capability)  than  the  Vtith  and  Clawson  [1972] 
curve  at  distances  less  than  about  40°,  because  at  these  distances  the  maximum  SNR  is  at  fre¬ 
quencies  greater  than  1  Hz.  We  simply  plot  the  Veith  and  Clawson  [1972]  curve  to  demon¬ 
strate  that  our  P-wave  detectability  estimates  are  generally  consistent  with  those  inferred  from 
the  Veith  and  Clawson  curve  at  longer  ranges,  where  the  dominant  signal  frequency  is  expected 
to  be  closer  to  1  Hz.  We  use  this  parameterization  to  represent  teleseismic  P-wave  detectabil¬ 
ity  for  all  network  c?”ibility  simulations  considered  in  this  report.  However,  we  recognize  that 
this  part  of  the  detectability  curve  is  more  uncertain  than  at  regional  distances,  so  we  do  not 
consider  networks  for  which  teleseismic  detection  capability  has  a  strong  influence  on  the  mag¬ 
nitude  threshold. 

We  were  not  able  to  parameterize  Lg  detectability  in  terms  of  the  frequency-dependent 
attenuation  arxl  noise  at  NORESS,  nor  have  we  separated  the  array  performance  from  the 
single-channel  capability.  Thus,  our  results  for  Lg  cannot  be  extrapolated  with  confidence. 
For  this  reason  we  only  include  Lg  in  simulations  for  networks  including  stations  assumed  to 
be  identical  to  NORESS.  That  is,  we  use  our  logarithmic  parameterization  (Figure  3.9)  to 
simulate  Lg  detectability  throughout  the  Soviet  Union,  though  we  include  Lg  blockage  by 
major  structural  boundaries. 

It  is  important  to  determine  the  sensitivity  of  our  detection  capability  estimates  to  our 
assumptions  about  Q.  Very  few  studies  of  regional  Pn  attenuation  have  been  conducted,  so  it 
is  difficult  to  estimate  the  amount  of  variability  from  region  to  region.  However,  numerous 
studies  of  Lg  attenuation  have  been  coivlucted.  Typically  the  1-Hz  Qi^  is  between  8(X)-11(X) 
in  the  eastern  United  States  arKl  the  exponent  of  a  power-law  frequency  dependence  is  between 
0.2  and  0.4  [e.g.,  Singh  and  Herrmann,  1983;  Hasegawa,  1985;  Gonez  et  al.,  1987;  Chun  et 
al,  1987;  Gupta  and  McLaughlin,  1987],  The  western  United  Sutes  is  characterized  by  a 
lower  1-Hz  Qt*  (140-200)  and  a  higher  hequetKy  exponent  (0.3-0.7)  [e.g.,  Singh  and 
Herrmann,  1983;  Nuttli,  1986;  Chavez  and  Priestley,  1986].  Figure  4.2  plots  predicted  Lg 
spectra  using  (2.1)  at  500  km  for  our  NORESS  Q  =  56(y®’*,  and  for  attenuation  models 
representative  of  the  eastern  and  western  United  States.  Also  plotted  are  predicted  spectra  for 
Q  =  84(y®'“  and  Q  =  28(y°'*,  which  correspond  to  ±50%  of  the  Qig  that  we  determined  for 
NORESS.  Note  that  this  range  in  Q  approximately  corresponds  to  the  difference  between 
spectra  recorded  in  the  eastern  and  western  United  States.  The  attenuation  we  observe  at 
NORESS  is  closer  to  attenuation  observed  in  the  eastern  United  States. 

Since  ±50%  of  our  NORESS  Q  seems  to  encompass  the  range  of  Qi^,  we  use  the  same 
range  of  models  for  Qp^  to  determine  the  sensitivity  of  our  network  simulations  to  that  parame¬ 
ter.  Figure  4.3  plots  predicted  Pn  spectra  at  5(X)  km  for  our  NORESS  Q  =  325/°  ^*,  and  for  Q 
equal  to  ±50%  of  that  value.  While  there  is  no  guarantee  that  this  range  is  representative  of 
the  variability  of  Pn  attenuation,  we  note  that  the  spectral  shapes  are  quite  dissimilar  and  imply 
very  different  conclusions  about  the  high-frequency  detection  capability. 

We  use  the  NORESS  noise  spectrum  for  all  rteiwoik  detection  capability  simulations 
presented  in  this  report.  It  is  trivial  to  determine  the  change  in  the  deteaion  capability  that 
would  result  from  a  shift  in  the  absolute  level  of  the  noise  spectrum.  Although  noise  estimates 
are  not  available  for  most  of  the  Soviet  Union,  Berger  et  al.  [198-8]  determined  average  rxrise 
levels  recorded  by  three  seismic  stations  in  the  vicinity  of  the  principal  imdetground  nuclear 
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Figure  4.2.  Predicted  Lg  spectra  for  a  magnitude  3.0  event  at  SOO  km  for  various  Q  models. 
The  solid  curves  are  for  Q(0  *  labeled  by  Go  The  spectnim  for  Go  =  560  is 

based  on  our  inverskm  of  NORESS  data  The  other  two  solid  curves  indicate  the  predicted 
spectra  if  G14  is  50%  greater  or  less  than  our  estimate  for  Scandiiuvia.  The  dashed  curves  are 
predicted  spectra  for  Q  estimates  typical  for  the  eastern  United  States  (EUS)  [Hasegawa,  1985] 
and  for  the  western  United  States  estimates  (WUS)  [Singh  and  Hermann,  1983]. 
Specifically,  we  use  G(/)  =  90Qr”  for  the  EUS  arai  G(/)  =  20(y°^  for  the  WUS.  The  bottom 
curve  is  the  average  ambient  NORESS  noise  spectrum. 
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test  site  in  East  Kazakh.  They  found  that  the  noise  levels  at  these  sites  are  roughly  equivalent 
to  those  observed  at  RSTN  (Regional  Seismic  Test  Network)  stations  between  1  and  20  Hz, 
and  10  to  20  dB  higher  than  observed  at  Lajitas.  Texas.  Suteau-Henson  and  Bache  [1988] 
found  that  NORESS  noise  levels  between  1  and  20  Hz  are  about  10-lS  dB  higher  than  those  at 
Lajitas.  Therefore,  the  NORESS  noise  levels  in  the  1-20  Hz  band  are  similar  to  those 
observed  in  East  Kazakh. 

It  is  also  of  interest  to  determine  the  extent  to  which  our  simulations  depend  on  the 
specific  characteristics  of  the  NORESS  signal  processing.  For  the  standard  beam  set  used  at 
NORESS,  we  find  that  the  array  offers  little  advantage  over  single  stations  at  high  frequency. 
However,  it  appears  that  the  NORESS  high-fiequency  beam  gain  could  be  improved  si^iply  by 
adding  mote  beams  [T.  Kvaema,  personal  communication,  1988].  Therefore,  we  simulate  the 
detection  capability  assuming  Vn  beam  gain  at  all  frequencies  for  all  stations,  and  compare  to 
simulations  normalized  by  the  observed  frequency-dependent  beam  gain  at  NORESS.  We  also 
simulate  the  detection  capability  of  a  network  composed  of  only  single  stations  (gain  =  1). 
These  two  extremes  bound  the  network  capability  for  any  combination  of  arrays  and  single  sta¬ 
tions.  We  also  note  that  some  improvement  in  detection  capability  may  be  possible  by  lower¬ 
ing  the  beam  thresholds,  albeit  a  cost  of  increased  false  alarms  to  be  recognized  by  the  pro¬ 
cessing  system  [Kvaerna  et  al.,  1987a]. 

Ftnally.  the  magnitudes  (A/J  used  in  our  regional  wave  detectability  study  are  based  on 
Lg  amplitude,  there  is  no  accurate  scaling  of  these  to  the  global  scale.  That  is,  we 
present  the  results  of  the  simulations  in  this  report  in  terms  of  NORESS  Af^  thresholds,  and 
these  may  not  represent  the  thresholds. 
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A2  Detection  capability  in  the  Soviet  Union 

This  section  presents  simulations  of  the  detection  capability  of  a  hypothetical  network  of 
seismic  stations/arrays  in  the  Soviet  Union  for  various  assumptions  about  the  characteristics  of 
regional  wave  propagation  in  the  USSR.  The  range  of  assumed  attenuation  models  and  station 
capability  is  based  on  our  experience  with  the  NORESS  array. 

An  event  detection  requires  detection  of  at  least  the  minimum  number  of  phases  neces- 
sary  to  locate  the  event.  For  regional  monitoring  using  arrays  and  three  component  instru¬ 
ments,  it  is  possible  to  estimate  both  azimuth  and  arrival  time  of  detected  phases.  The  determi¬ 
nation  of  epicenter  and  origin  time  require  arrival  time  aiul  azimuth  information  from  at  least 
two  phases,  Depth  estimation  is  pooriy  constrained  by  arrival  time  data  for  crustal  events  at 
regional  distances.  Well-constrained  depth  estimates  require  the  detection  of  depth  phases,  but 
this  is  difScult  to  represent  in  the  simulations.  Thus,  in  all  examples  to  be  presented,  we 
require  that  at  least  three  phases  (arrival  times  and  azimuths)  are  detected  by  the  hypothetical 
network. 

We  use  a  network  consisting  of  13  stations  external  to  the  Soviet  Union  (sited  at  the 
location  of  existing  stations)  and  20  equally  spaced  stations  internal  to  the  Soviet  Union  for  all 
of  our  detection  capability  simulations  (Figure  4.4).  The  inter-station  spacing  for  the  internal 
network  is  approximately  1000  km.  Braa  et  al.  [1987]  show  that  for  this  network,  the  location 
uncertainty  for  events  near  the  network  threshold  for  3  detections  is  about  10-15  km  (this 
uncertainty  can  often  be  reduced  with  master  event  location  techruques). 

With  this  many  internal  stations,  the  detection  capability  is  only  weakly  dependent  upon 
our  estimates  of  teleseismic  P-wave  detectability.  For  example,  simulations  that  eliminate  all 
detections  beyond  12°  produce  90%  thresholds  that  differ  by  less  than  0.1  from  those 
obtained  when  detections  are  allowed  to  teleseismic  ranges. 

Table  4.2  lists  the  parameters  and  summarizes  the  results  of  the  simulations  presented  in 
this  report.  Simulations  1-4  are  normalized  by  our  temporal  parameterization  of  detectability 
(Figures  3.7  and  3.9).  These  simulations  assume  that  the  signal  and  noise  characteristics  in  the 
USSR  are  identical  to  those  at  NORESS,  and  that  each  station  is  a  NORESS-type  array  with 
the  standard  beam  set  listed  in  Table  3.1.  Figure  4.S  plots  the  50%  and  90%  Ml  thresholds  for 
detecting  three  Pn  phases.  The  90%  Ml  threshold  is  between  2.4  and  2.7  for  events  in  the 
Soviet  Union.  The  50%  dueshold  is  between  2.2  and  2.5.  Figure  4.6  plots  the  50%  and  90% 
Ml  thresholds  for  detecting  three  P  or  Lg  phases  for  the  same  netwoik.  The  90%  Ml  threshold 
is  between  2.1  and  2.6  for  events  in  the  Soviet  Union.  The  50%  threshold  is  between  1.8  and 
2.3.  Since  we  cannot  extrapolate  Lg  detectability  to  regions  with  different  attenuation,  we  do 
not  include  Lg  in  the  detection  criteria  for  the  remaining  simulations.  We  simply  note  that 
based  on  NORESS  observations,  we  expect  the  Ml  threshold  to  be  0.2  to  0.3  lower  for  this 
network  if  Lg  is  included. 
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SIMULATION  NETWORK 
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Figure  4.4.  Location  of  seismic  sutions  comprising  the  hypotheUcal  networic  used  in  our 
detection  capabiUty  simulaUons.  Open  circles  denote  sites  to  which  Lg  propagation  from  the 
USSR  is  prohibited. 
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Figure  4.5.  Contours  of  the  (a)  50%  and  (b)  90%  Af/,  threshold  for  detecting  3  Pn  phases  for  a 
nctwoik  of  NORESS-type  arrays;  13  external  to  the  Soviet  Union  ani  20  intenud.  It  is 
assumed  that  attenuation,  noise,  and  array  perfonnance  are  all  identical  to  that  observed  at 
NORESS.  The  nonnalizadon  is  based  on  our  temporal  parameterization  of  Fn  detectabiUty. 
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Figure  4.6.  Coniouis  of  the  (a)  50%  and  (b)  90%  threshold  for  the  detection  of  3  Pn  or  Lg 
phases  for  a  network  of  NORESS-type  arrays;  13  external  to  the  Soviet  Union  and  20  internal 
It  is  assumed  that  attenuation,  noise,  and  array  performance  are  all  identical  to  that  observed  at 
NORESS.  The  normalization  is  based  on  our  temporal  parameterization  of  Pn  and  Lg  detecu- 


Table  4.2.  Simulations 


Simulation 

Number 

Detection 

Criterion 

Gain 

Qpn 

Qo 

n 

Parameterization 

Ml 

Threshold 

1 

3  F-90% 

NORESS 

Temporal 

2 

3  /»-50% 

NORESS 

- 

— 

Temporal 

2.2-2.5 

3 

3  (P  or  Lg)-90% 

NORESS 

— 

- 

Temporal 

2. 1-2.6 

4 

3  iP  or  Lgy50% 

NORESS 

- 

- 

Temporal 

1. 8-2.3 

5 

3P-90% 

NORESS 

325 

0.48 

Spectral 

2.4-2.7 

6 

3  P-m% 

Vn 

325 

0.48 

Spectral 

i-y-in 

7 

3  /*-90% 

1.0 

325 

0.48 

Spectra] 

2.7-3.3 

8 

3  P-90% 

NORESS 

487.5 

0.48 

Spectral 

2.3-2.6 

9 

3  P-90% 

NORESS 

162.5 

0.48 

Spectral 

2.1-29 

Simulations  S-9  of  Table  4.2  use  our  spectral  parameterization  of  Pn  detectability.  Figure 
4.7  shows  the  corresponding  Pn  detectability  curves.  Note  that  these  curves  include  an  imi^cit 
frequency  dependence,  since  at  each  distance  they  represent  the  detecubility  of  the  beam  with 
the  maximum  SNR.  This  beam  is  determined  by  a  combination  of  the  single-channel  SNR  and 
the  hequency-dependent  beam  gain.  For  example,  our  spectral  model  predicts  the  frequency  of 
the  maximum  SNR  on  a  single-channel  to  be  greater  than  8  Hz  for  distances  less  than  400-500 
km.  However,  the  optimal  detectability  is  predicted  to  be  on  the  3-S  Hz  beam  since  the  stan¬ 
dard  NORESS  beam  set  gives  gains  close  to  one  at  high  frequency. 

Simulation  S  uses  our  NORESS  Pn  attenuation  model  and  the  beam  gains  observed  at 
NORESS  (e.g.,  Vn  between  2-S  Hz.  and  one  for  the  4-8  Hz  and  8-16  Hz  beams).  The  90% 
Ml  threshold  for  detecting  three  P  phases  is  shown  in  Figure  4.8.  We  use  the  attenuation  and 
source  parameters  derived  from  the  inversion  of  NORESS  spectra,  the  average  ambient 
NORESS  noise  spectrum,  and  the  beam  gains  observed  at  NORESS.  Note  that  the  result  is 
nearly  the  same  as  the  simulation  normalized  by  our  temporal  parameterization  of  NORESS  Pn 
detectal^ty.  This  demonstrates  the  internal  consistency  of  our  temporal  and  ^rectral  parame- 
terizations. 

Figure  4.9a  plots  the  90%  M^  thresholds  for  detecting  three  P  phases  for  a  network  of 
arrays  with  Vn  beam  gain  at  all  frequetKies.  Attenuation  and  source  parameters  derived  from 
NORESS  spectra  and  the  average  NORESS  noise  spectrum  were  used  to  calculate  detectability. 
For  the  network  of  arrays,  we  assume  the  same  sub-arrays  are  used  in  beamforming  as  those 
listed  for  the  standard  NORESS  beams  in  Table  3.1  (e.g.,  nine  elements  are  used  for  the  two 
highest-frequency  beams  and  up  to  22  elements  are  used  for  the  lower  frequency  beams).  We 
find  diat  the  optimal  detectability  is  on  the  8-16  Hz  beam  to  a  distance  of  3°  and  on  the  4-8  Hz 
beam  from  3®  to  12®.  The  90%  M^  threshold  is  between  2.3  and  2.7  for  events  in  the  Soviet 
Union.  Therefore,  by  optimizing  beam  gain  at  high-frequency,  it  is  possible  to  reduce  the  M^ 
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Figure  4.7.  Predicted  regional  Pn  detectability  curves  based  on  (3.4)  and  (3.5)  with  a  =  1.35 
for  several  assumptions  about  Q  and  beam  gaia  The  curves  with  Qf,^  s  325/°^^  illustrate  the 
dependence  of  Pn  detecubility  on  beam  gain  (the  beam  gain  observed  at  NORESS  for  die 
standard  beam  set  is  Vn  between  2-5  Hz  and  one  for  the  two  hi^r  frequency  beams).  The 
dependence  of  Pn  detectability  on  Q  is  demonstrated  by  the  curves  with  closed  symbols. 
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Rgure  4.9.  Contours  of  the  90%  direshold  for  detecting  3  Pn  phases  for  a  netwoik  of  (a) 
NORESS-type  arrays  with  Vyv  beam  gain  at  ail  frequencies  and  (b)  single-stations.  It  is 
assumed  that  attenuation  and  noise  are  identical  to  that  observed  at  NORESS.  The  normaliza¬ 
tion  is  based  on  our  spectral  parameterization  of  Pn  detectability. 


detection  threshold  by  about  0.1  throughout  most  of  the  Soviet  Union.  For  comparison.  Figure 
4.9b  plots  the  90%  Mi  threshold  for  a  network  composed  of  single  stations.  To  simulate  the 
detection  c^)ability  of  single  stations,  we  simply  set  the  beam  gain  to  one  for  all  frequencies. 
The  Pn  detectability  is  lower  than  that  for  the  array  network  by  about  O.S  at  all  distances  (Fig¬ 
ure  4.7).  However,  since  log  moment  increases  at  a  rate  of  l.SSAf^,  the  difference  in  the  Mi 
threshold  is  slightly  less  than  this.  The  90%  Mi  threshold  for  the  network  of  single  stations  is 
about  2.7  to  3.3  for  events  in  the  Soviet  Union. 

Our  final  examples  are  simulations  of  the  detection  capability  in  the  Soviet  Union  of  a 
network  of  NORESS-type  arrays  with  50%  higher  and  lower  than  we  estimate  Scandina¬ 
via.  The  90%  Ml  thresholds  for  detecting  three  Pn  phases  are  shown  in  Figure  4.10.  We  use 
the  average  NORESS  noise  spectrum  and  the  beam  gains  observed  for  the  standard  NORESS 
beam  set  The  90%  Mi  threshold  if  Q  is  50%  higher  than  at  NORESS  is  between  2.3  to  2.6 
for  events  in  the  Soviet  Union.  If  Q  is  50%  lower,  the  threshold  is  between  2.7  and  2.9.  That 
is,  the  threshold  changes  by  about  0.4  over  diis  range  of  Q. 

The  accuracy  of  our  simulations  of  the  detection  capability  of  hypothetical  networks  in 
the  Soviet  Union  is  governed  by  the  accuracy  of  our  estimates  of  attenuation  and  noise  in 
region.  We  are  confident  that  we  have  accurately  tepmented  the  detection  capabUiO'  of 
NORESS  array,  however  it  is  not  possible  to  assess  the  accuracy  of  our  estimates  for 
Soviet  Union  without  examining  data  collected  in  that  region. 
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Figure  4.10.  Contours  of  the  90%  threshold  for  detecting  3  Pn  phases  for  a  nttwoik  of 
NORESS-type  arrays  assuming  that  Qp^  is  (a)  50%  higher  and  (b)  50%  lower  than  our 
NORESS  estimate.  It  is  assumed  that  array  performance  is  identical  to  that  observed  at 
NORESS.  The  normalization  is  based  on  our  spectral  parameterization  of  Pn  detectability. 


5.  CONCLUSIONS 


The  main  objective  of  this  study  is  to  evaluate  detection  capability  of  a  hypothetical  net- 
woik  of  seismic  stations/arrays  in  and  around  the  Soviet  Union.  We  use  a  parameterization  of 
regional  wave  spectra  recorded  at  NORESS  as  a  basis  for  nonnalizing  simulations  of  network 
capability.  Since  accurate  attenuation  and  noise  estimates  are  not  available  for  the  Soviet 
Union,  we  also  detennine  the  sensitivity  of  our  results  to  our  assumptions  about  the 
frequency-dependent  signal  and  noise  characteristics  in  the  USSR. 

Our  approach  involves  four  steps.  First,  we  parameterize  NORESS  spectra  in  terms  of  Q 
and  seismic  moment.  Next,  we  determine  the  relationship  between  the  sirtgle-chanrKl  spectral 
amplitudes  and  the  temporal  amfditudes  used  in  sigiud  detection  at  NORESS.  Third,  we 
predict  the  detection  crqMtbili^  of  the  NORESS  array  based  on  our  spectral  parameterization 
and  validate  our  results  by  comparing  to  observed  detection  statistics  and  to  the  results  of  an 
empirical  study  of  detection  capability  in  this  regkm.  Firudly,  we  extrapolate  our  results  for 
NORESS  to  simulate  the  detection  capability  of  a  hypothetical  network  of  seismic  stations  in 
and  around  the  Soviet  Union. 

From  our  parameterization  of  NORESS  spectra,  we  conclude: 

•  The  attenuation  of  Pn  spectra  at  NORESS  between  I-IS  Hz  is  consistent  with  r~^-^ 
geometric  spreading  and  Q(f)  = 

•  The  attenuation  of  Lg  spectra  at  NORESS  between  1-7  Hz  (computed  for  the  fixed 
group  velocity  window,  3.6-3.0  km^)  is  consistertt  with  cylmdrical  spreading 
beyond  a  transition  distance  of  KX)  km,  and  Q{f)  = 

These  results  were  obtained  by  inverting  spectra  for  186  regiorutl  events  with  b^een  l.l 
and  4.8.  The  inversion  produced  a  range  of  attenuation  models  consistent  with  the  observa¬ 
tions.  To  resolve  parameter  trade-offs  we  impose  the  constraint  that  the  separate  Pn  and  Lg 
inversions  define  an  internally  consistent  set  of  smrroe  parameters.  This  is  the  basis  for  select¬ 
ing  our  "preferred"  model.  However,  all  the  models  give  similar  predictions  for  detection 
capability  for  the  distance  and  frequency  range  of  the  NORESS  data.  Our  model  of  Pn 
attenuation  at  NORESS  imidies  the  following  conclusion  about  network  detection  capability: 

•  If  the  upper  mantle  structure  of  the  Russian  platform  is  at  all  similar  to  Scandinavia, 
high  frequencies  (>20  Hz)  are  expected  to  contribute  very  little  to  network  detection 
capalMlity,  except  at  distances  less  than  about  4(X)-S00  km.  At  these  distances,  the 
detection  thresholds  are  quite  low  in  lower  frequency  bands  (about  2.1  for  Pn  detec¬ 
tion  and  1.8  for  detection  of  either  Pn  or  Lg  at  NORESS).  Therefore,  our  results 
for  NORESS  are  inconsistent  with  the  major  conclusion  of  Evernden  at  al.  [1986]; 
that  detection  capability  for  events  in  the  USSR  would  be  enhanced  at  hequencies 
greater  than  20  Hz. 
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To  use  our  parameterization  of  regional  Pn  spectra  to  predict  and  extrapolate  NORESS 
detection  capability,  we  determine  the  relationship  between  the  single-channel  spectral  ampli¬ 
tudes  and  the  time  domain  amplitudes  used  in  signal  detection.  Our  main  results  are: 

•  The  Pn  SNR  used  in  signal  detection  at  NORESS  can  be  expressed  as  a  product  of 
the  single-channel  spectral  SNR  and  terms  specific  to  the  NORESS  array 
configuration  and  beamforming.  This  factorization  allows  us  to  predict  Pn  detecta¬ 
bility  for  other  regions,  or  for  other  station  configurations. 

•  Due  to  dispersion  and  the  non-stationarity  of  the  pte-Lg  noise,  a  similar  relationship 
could  not  be  found  for  Lg.  Therefore,  we  simply  determir^  an  empirical  relation¬ 
ship  for  the  Lg  detectability  at  NORESS  and  note  that  our  results  cannot  be  extrapo¬ 
lated  with  confidence  to  regions  with  different  attenuation. 

From  our  parameterization  of  Pn  spectra,  and  the  relationship  between  temporal  and  spec¬ 
tral  amplitudes,  we  estimate  the  detection  capability  of  the  NORESS  array.  We  also  estimate 
Pn  and  Lg  detectability  at  NORESS  based  on  a  simple  parameterization  of  observed  temporal 
amplitudes.  We  conclude  the  following  for  the  detection  capability  at  NORESS: 

•  The  S0%  and  90%  thresholds  for  detecting  Pn,  based  on  our  parameterization  of 
NORESS  spectra,  are  1.6  and  2.0  at  400  km  and  2.5  and  2.8  at  1000  km. 

•  Based  on  our  parameterization  of  observed  temporal  amplitudes,  the  50%  and  90% 
Ml  thresholds  for  detecting  Pn  are  1.6  and  2.0  at  400  km  and  2.4  and  2.8  at  1000 
km.  These  are  consistent  with  tlrc  estimates  obtained  from  our  spectral  parameteri¬ 
zation.  Our  estimates  of  the  50%  and  90%  Mi  thresholds  for  detecting  either  Pn  or 
Lg  at  1000  km  ate  2.2  and  2.6.  Thus,  including  Lg  in  the  detection  criterion 
reduces  the  magnitude  threshold  by  about  0.2  at  1000  km. 

•  Our  results  for  NORESS  detection  capability  are  consistent  with  those  of  Ringdal 
[1986],  who  compared  detections  at  NORESS  with  bulletins  produced  by  local 
seismic  networks.  His  estimates  of  the  50%  and  90%  Mi  thresholds  between  700- 
1400  km  are  2.3  and  2.7  for  Pn  detection  and  1.9  and  2.5  for  either  Pn  or  secon¬ 
dary  phase  detectioa 

We  are  confident  that  we  have  accurately  estimated  the  detection  capability  of  the  NORESS 
array.  However,  these  estimates  ate  based  on  the  standard  beam  set  used  to  process  NORESS 
data.  For  this  set,  we  found  that  close  to  Viv  beam  gain  (where  N  is  the  number  of  elements 
used  in  beamforming)  is  achieved  between  2-5  Hz,  but  that  the  gain  at  high  frequency  (8-16 
Hz)  is  close  to  one.  It  is  possible  that  the  beam  gain  could  be  improved  at  high  frequency  by 
simply  adding  steered  beams  [T.  Kvaema,  personal  communication,  1988].  If  this  is  the  case, 
our  results  may  underestimate  the  detection  capability  of  the  NORESS  array  with  an  optimal 
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beam  dq)loymeiU.  Of  course,  this  is  only  relevant  for  distances  less  than  about  500  km.  since 
at  longer  ranges  the  single-channel  SNR  is  maximum  at  lower  hnequency. 

The  final  results  of  this  study  are  simulations  of  the  detection  capability  of  a  hypothetical 
network  in  the  Soviet  Union.  The  network  consists  of  13  stations  external  to  the  USSR  and  20 
internal  stations.  The  spacing  of  the  internal  stations  is  about  1000  km.  We  did  not  simulate 
the  capability  of  sparse  networks  due  to  greater  uncertaituy  in  the  normalization  at  teleseismic 
distances.  We  conclude  the  following: 

•  First,  we  assume  that  the  frequency-dependent  signal  and  noise  characteristics  in  the 
Soviet  Union  are  identical  to  those  we  observe  at  NORESS.  and  that  the  monitoring 
network  consists  of  NORESS-type  arrays  equipped  with  the  standard  beam  set  used 
to  process  NORESS  data  (beam  gain  equal  to  '<In  between  2-S  Hz.  and  one  at 
higher  frequency).  Based  on  these  assumptions,  the  90%  threshold  for  detecting 
3  Pn  phases  for  events  in  the  Soviet  Union  is  between  2.4  and  2.7.  Including  Lg  in 
the  detection  criterion  reduces  the  foreshold  by  0.2-0.3. 

•  To  determine  the  extent  to  which  our  results  depend  on  the  array  performatKe 
observed  at  NORESS.  we  assume  that  the  frequency-dependent  signal  and  raise 
characteristics  in  the  Soviet  Union  are  identical  to  those  at  NORESS.  but  we  vary 
assumptions  about  beam  gain.  The  90%  direshold  for  detecting  3  Pn  phases  for 
a  network  of  NORESS-type  arrays  with  beam  gain  at  all  frequencies  (i.e.,  detec¬ 
tion  of  high  frequencies  better  than  that  currently  observed  at  NORESS)  is  between 
2.3-2.1.  Since  the  enhanced  SNR  gain  at  high  frequencies  is  relevant  only  for  dis¬ 
tances  less  than  about  500  km.  the  threshold  for  events  near  the  periphery  of  the 
internal  network  is  insensitive  to  the  high-frequency  beam  gaia 

•  The  Ml  threshold  for  a  network  of  single  stations  is  between  0.4  and  0.5  higter 
than  that  for  the  network  of  NORESS-type  arrays  with  beam  gain  at  all  frequen¬ 
cies.  It  is  0.3  to  0.5  higher  than  for  a  network  of  arrays  with  observed  NORESS 
capabilities. 

•  If  Qp^  is  50%  higher  in  the  Soviet  Union  than  our  estimate  for  Scandinavia,  the 
90%  Ml  threshold  for  detecting  3  Pn  phases  for  a  i^twork  of  NORESS-type  arrays 
is  between  2.3  and  2.6.  The  90%  Ml  threshold  is  between  2.7  and  2.9  if  Qp„  is 
50%  lower  than  our  estimate  for  Scandinavia.  Very  few  studies  of  frequency- 
dependent  Pn  attenuation  have  been  conducted,  so  it  is  difficult  to  bound  the 
expected  regional  vaiiaUlity  of  Qp^  The  ±50%  range  was  based  on  comparing  our 
Qi^  to  values  estimated  for  the  eastern  and  western  United  States.  While  there  is 
no  guarantee  that  this  range  is  representative  of  the  variability  of  Pn  attenuation,  the 
spectral  shapes  are  quite  dissimilar  and  result  in  very  different  exclusions  about 
high-frequency  detection  capability. 
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The  accuracy  of  these  simulations  of  the  detection  capability  of  hypothetical  networks  depends 
on  the  accuracy  of  our  estimates  of  attenuation  and  noise  in  the  Soviet  Union.  Extrapolation  of 
observations  from  other  regions  (such  as  the  region  around  NORESS  used  in  this  study)  have 
large  uncertainty.  To  be  more  confident,  we  need  to  study  data  from  stations  throughout  the 
Soviet  Union  to  estimate  the  actual  attenuation  and  noise  conditions  to  be  encountered. 
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